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Preface 


The  mathematical  model  described  In  this  paper  Is 
Intended  primarily  for  use  by  systems  engineers,  and  for 
that  reason  I  have  adopted  to  a  degree  the  viewpoint  of  the 
systems  engineer  rather  than  that  of  the  theoretical 
physicist  In  Its  development.  I  have  made  a  determined 
effort  to  describe  the  four-level  Laser  system  In  terms  of 
parameters  which  are  physically  measvirable  In  the  laboratory 
by  means  of  relatively  simple  and  stral^tforward  methods. 
Whether  or  not  such  an  approach  Is  adequate  to  the  problem 
at  hand  remains  to  be  detezmlned  through  the  actual 
application  of  the  model  to  physical  systems. 

I  wish  to  express  my  gratitude  to  the  many  Individuals 
who  have  assisted  me  In  this  study.  In  particular,  I  wish 
to  thank  Mr.  Slrons,  of  the  Navigation  and  Guidance 
Laboratory,  Aeronautical  Systems  Division,  who  sponsored  the 
project  and  supplied  the  equipment  for  the  study  of  xenon 
lamp  pulses  and  the  sample  for  spectroscopic  analysis,  Mr. 
Blaslngame,  of  the  Electronic  Technology  Laboratory,  ASL, 
who  performed  the  spectroscopy,  and  Professor  Lubelfeld,  of 
the  Electrical  Engineering  Department,  Air  Force  Institute 
of  Technology,  my  Faculty  Thesis  Adviser. 

Tom  G.  Purnhagen 
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Abstract 


The  Air  Force  Is  considering  several  pulsed  four-level 
Lasers  for  possible  use  in  weapons  systems.  So  far,  the 
advantages  of  these  Lasers  over  the  ruby  Laser,  such  as 
lower  threshold  and  reduced  spiking,  have  been  offset  by 
their  lower  efficiency.  Although  four-level  Laser  efficiency 
could  be  maximized  by  experimental  means,  analysis  by  means 
of  a  mathematical  model  appears  more  suitable  for  economic 
reasons.  Such  a  model  is  developed  in  this  paper  by 
identifying  the  loss  mechanisms  >diich  operate  in  the  four- 
level  Laser  system,  defining  the  system  efficiency  as  the 
product  of  six  sub-efficiencies,  one  associated  with  each 
loss  mechanism,  and  deriving  equations  for  the  sub- 
efficiencies  in  terms  of  physically  measurable  parameters 
of  the  Laser  system.  The  model  is  limited  in  its  present 
form  to  the  representation  of  pulsed  four-level  Lasers 
employing  helical  or  annular  lamp  pumping,  but  could  be 
easily  extended  to  cover  other  pumping  geometries  and 
continuous  operation. 

Based  on  the  mathematical  model,  a  program  is  written 
for  the  IBM  1620  digital  computer  ^ich  calculates  the 
efficiency  of  a  pulsed  four-level  Laser  oscillator  or  the 
maximum  effiolenoy  and  maximum  gain  of  a  pulsed  four-level 
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Laser  amplifier.  The  program  Is  written  In  a  modification 
of  IBM  FORTRAN,  and  consists  of  three  sub-programs;  one 
^Ich  need  only  be  run  once,  one  vdilch  must  be  run  once  for 
each  Laser  material,  and  one  which  produces  the  final 
results.  Output  data  from  part  1  are  listed,  and  detailed 
operating  Instructions  for  parts  2  and  3  are  given.  A 
sample  problem  Involving  the  design  of  a  neodymium  doped 
glass  Laser  oscillator  Is  solved  to  Illustrate  the  use  of 
the  model.  On  the  basis  of  the  model  and  the  experience 
gained  In  Its  develcpment  and  testing.  It  can  be  concluded 
that  the  efficiency  of  pulsed  four-level  Lasers  will  tend 
to  Increase  as  length,  radius,  doping  density,  and  pumping 
energy  are  Increased,  and  decrease  as  end  reflectivities 
and  pumping  pulse  time  constant  are  Increased,  within  the 
limits  Imposed  by  the  assumptions  and  approximations  used 
In  the  development  of  the  model. 

The  results  produced  by  the  model  have  not  been  checked 
experimentally,  and  It  Is  recommended  that  such  a  check  be 
made  before  the  model  Is  used  In  actual  design  problems. 
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A  OOlCPUTER  MODEL  FOR  THE  DETERMINATION 
OF  FOUR-LEVEL  LASER  EFFICIEirOY 


I.  Introduction 


The  United  States  Air  Force  Is  presently  conducting 
research  to  determine  the  potentialities  of  the  Laser  as  a 
military  weapon.  Although  much  of  the  current  effort  Is 
centered  about  the  ruby  Laser,  which  ifas  the  first  solid 
state  Laser  developed  and  hence  Is  the  most  highly 
perfected,  the  Interest  of  the  Air  Force  Is  by  no  means 
limited  to  that  particular  device.  Numerous  other  materials 
are  also  being  Investigated,  among  >dilch  are  various 
crystals,  glasses,  and  liquids  doped  with  divalent  or 
trlvalent  rare  earth  Ions. 

The  rare  earth  Lasers  operate  In  a  four-level  energy 
configuration  (Figure  1b)  as  opposed  to  the  three-level 
ruby  Laser  (Figure  la).  A  comparison  of  the  two  energy 
level  configurations  Indicates  that  the  four-level  system 
has  some  Inherent  advantages  over  the  three-level  system, 
among  >dilch  are  a  lower  threshold  and  reduced  spiking 
(relaxation  oscillations).  For  high  power  applications, 
however,  these  advantages  have  been  offset  by  the  relative 
Inefficiency  of  existing  four-level  systems  compared  with 
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(a)  Three-level  System,  (b)  Pour-level  System. 

the  ruby.  If  the  rare  earth  Lasers  are  to  be  useful  in 
weapons  systems,  then,  It  Is  Important  to  be  able  to  design 
and  operate  them  in  the  most  efficient  manner  possible. 

The  problem  of  maximizing  Laser  efficiency  can  be  (and 
has  been)  approached  through  experimental  study;  that  is, 
through  trying  various  combinations  of  physical  parameters, 
doping  levels,  end  mirror  reflectivities,  and  so  forth,  to 
see  yrhioh  combination  gives  the  best  results.  This  method, 
although  practical,  is  expensive  and  time  consuming.  A 
better  approach  would  seem  to  be  to  develop  a  mathematical 
model  for  the  Laser  system,  and  to  analytically  (or 
graphically)  determine  the  optimum  design  prior  to  the 
actual  physical  construction  of  the  Laser. 

The  remainder  of  this  paper  is  devoted  to  the 
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development  of  a  program  for  the  IBM  1620  digital  computer 
iidiloh  calculates  the  efficiency  of  a  pulsed  four-level  Laser 
from  the  values  of  the  system  parameters  and  certain 
physically  measurable  characteristics  of  the  Laser  material 
and  the  pumping  source.  Although  the  program  and  the 
mathematical  model  from  which  It  Is  derived  were  originally 
designed  to  represent  Lasers  employing  trlvalent  neodymium, 
It  Is  believed  that  the  results  are  applicable  directly  to  a 
broad  class  of  four-level  Lasers,  and  that  only  minor 
changes  would  be  required  to  further  extend  the  range  of 
applicability. 

Chapter  II  of  this  paper  Identifies  the  loss  mechanisms 
In  the  four-level  Laser  system  and  defines  the  overall 
efficiency  In  terms  of  six  sub-efficiencies.  Chapter  III 
lists  the  more  Important  assumptions  and  approximations  on 
which  the  analysis  Is  based,  and  Indicates  how  these  will 
affect  the  range  of  applicability  of  the  model.  The  actual 
development  of  the  mathematical  model  Is  presented  In 
Chapter  IV,  and  In  Chapter  V  this  model  is  translated  Into 
the  computer  program.  Chapter  VI  Illustrates  the 
application  of  the  computer  program  to  a  design  problem, 
and  Chapter  VII  Is  a  summary  of  the  results  and  conclusions 
irtxloh  may  be  drawn  from  this  study. 


3 


GE/EE/62-18 


II.  The  Loss  Meohaniems 

Before  a  model  can  be  developed  to  determine  the 
efficiency  of  the  four-level  lAser  system,  it  is  necessary 
to  identify  the  loss  mechanisms  which  are  operative  in  the 
system.  This  can  be  done  by  following  the  energy  through 
the  system  from  input  to  output  and  noting  the  losses  which 
occur. 

Lamp  Losses 

Nearly  all  the  Laser  systems  now  in  use  employ  some 
sort  of  transducer  which  converts  electrical  energy  into 
photon  energy  to  excite  the  active  ions  in  the  Laser  itself. 
In  all  such  transducers  there  is  some  loss  of  energy  due  to 
the  inefficiency  of  the  conversion  process.  Accordingly,  a 
quantity  the  lamp  efficiency,  can  be  defined  as  the 
ratio  of  the  output  energy  of  a  given  pumping  source  which 
lies  within  a  specified  range  of  wavelength  to  the  total 
input  energy  to  that  source. 

Geometrical  Losses 

The  portion  of  the  photon  energy  produced  by  the 
pumping  source  which  actually  enters  the  Laser  rod  will 
depend  on  the  geometry  of  the  lamp  and  enclosure,  and  to 
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some  degree  on  other  factors  such  as  the  refractive  Index  of 
the  Laser  material  and  the  reflectivity  of  the  walls  of  the 
enclosure.  The  remainder  of  the  lamp  output  Is  lost  energy. 
For  any  particular  pumping  scheme ,  a  quantity  the 
geometrical  efficiency,  can  he  defined  as  the  ratio  of  the 
photon  energy  which  enters  the  Laser  rod  to  that  produced  by 
the  pumping  source,  again  within  a  specified  wavelength 
range . 

Capture  Losses 

Part  of  the  energy  which  enters  the  Laser  rod  Is  used 
to  pump  Ions  from  the  ground  state  to  the  pumping  band.  The 
rest  of  the  energy  passes  through  the  rod  and  Is  lost.  The 
amount  of  energy  lost  In  this  manner  depends  on  the 
dimensions  of  the  rod,  the  density  of  active  Ions,  and  the 
degree  to  which  the  absorption  spectnun  of  the  active 
material  matches  the  spectrum  of  the  pumping  source.  The 
ratio  of  the  energy  transferred  to  active  Ions  to  the  total 
energy  entering  the  rod  will  be  defined  as  the  capture 
efficiency,  Aq,  of  the  system. 

Relaxation  Losses 

When  the  active  Ions  decay  from  the  pumping  band,  some 
of  their  energy  Is  expended  In  the  desired  transition  from 
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the  metastable  state  to  the  terminal  state,  idille  the  rest 
of  the  energy  Is  dissipated  In  non-radlatlve  transitions  or 
radiated  at  frequenoles  other  than  the  Laser  output 
frequency.  Ihe  ratio  of  the  energy  radiated  In  the  desired 
transition  to  the  total  energy  stored  In  the  pumping  band 
will  be  'defined  as  the  relaxation  efficiency,  y].Q.  This 
quantity  Is  equal  to  the  quantum  efficiency  £  (the  ratio  of 
the  number  of  photons  radiated  at  the  desired  frequency  to 
the  total  number  of  photons  absorbed  In  the  pumping  process) 
multiplied  by  the  ratio  of  the  energy  of  an  output  photon  to 
the  average  energy  of  the  pumping  photons. 

Radiation  Losses 

^art  of  the  energy  available  In  the  desired  transition 
goes  Into  the  coherent  radiation  field,  and  the  rest  Is  lost 
through  spontaneous  emission.  The  magnitude  of  the 
radiation  loss  depends  on  how  quickly  threshold  Is  reached 
and  on  the  fluorescent  lifetime  of  the  metastable  state. 

The  corresponding  radiation  efficiency,  Hr*  is  defined  as 
the  ratio  of  the  energy  which  enters  the  coherent  field  to 
the  total  energy  involved  in  the  transition  from  the 
metastable  state  to  the  terminal  state. 
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Optical  Losses 

Finally,  some  of  the  energy  In  the  coherent  field  Is 
channeled  Into  useful  output,  \dalle  the  rest  Is  absorbed  or 
scattered  by  the  Laser  rod  and  the  end  mirrors.  The  optical 
efficiency,  IQq,  can  be  defined  as  the  ratio  of  useful  output 
energy  to  the  total  energy  idilch  enters  the  coherent  field. 

Overall  Efficiency 

If  the  lamp  efficiency,  geometrical  efficiency,  capture 
efficiency,  relaxation  efficiency,  radiation  efficiency,  and 
optical  efficiency  of  a  Laser  system  are  defined  as  stated 
In  the  preceding  paragraphs,  then  the  overall  system 
efficiency  Is  simply  the  product  of  the  six  sub-efflclenoles, 

(2-1) 
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III*  Ab sumptions  and  Approximations 

In  order  to  keep  the  model  from  beoomlng  enormously 
complicated,  certain  assumptions  and  approximations  will  be 
made  concerning  the  system  and  Its  representation  by  the 
model.  These  assvmptlons  and  approximations  fall  Into 
three  general  categories:  (1)  Assumptions  about  the  geometry 
and  physical  characteristics  of  the  system,  (2)  assumptions 
concerning  the  mathematical  foxmxilatlon  of  the  efficiency 
equations,  and  (3)  numerical  approximations  In  the  computer 
program. 

Assumptions  About  the  System 

The  Laser  Rod.  The  active  element  In  the  Laser  system 
will  be  assumed  to  be  a  right  olroular  cylinder  of  some 
homogeneous  material  doped  uniformly  over  Its  volume  with 
€Ln  appropriate  active  Ion.  The  ends  of  the  rod  will  be 
assumed  to  be  perfectly  flat  euid  perpendicular  to  the  axis 
of  the  cylinder. 

Illvunlnatlon.  It  will  be  assumed  that  the  Laser  rod 
Is  Illuminated  imlformly  over  the  cylindrical  surface,  that 
the  effects  of  reflection  and  refraction  at  the  surface  can 
be  neglected,  and  that  the  ll£^t  flux  entering  the  rod 
through  any  element  of  surface  area  dA  Is  diffused  evenly 
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over  a  solid  angle  of  2it  steradlans  within  the  rod.  It  Is 
further  assumed  that  the  flux  density  In  the  rod  Is 
Independent  of  distance  along  the  axis  and  that  the  flux 
density  can  be  calculated  as  If  the  rod  were  Infinitely 
long,  Implying  that  either  the  rod  length  Is  much  greater 
than  the  rod  radius,  or  both  ends  have  high  reflectivity,  or 
both. 

It  Is  believed  that  the  above  assumptions  are 
reasonably  valid  for  most  solid  state  Lasers  pumped  by 
helical  or  annular  lamps.  The  use  of  the  model  to  represent 
other  pimping  schemes  may  require  some  modification  of  the 
first  part  of  the  computer  program  to  adjust  for  the 
difference  In  flux  distribution  In  the  rod. 

Assumptions  Used  In  the  Mathematical  Development 

The  model  described  In  this  paper  Is  Intended  to  be 
applicable  to  a  wide  variety  of  Laser  systems,  and 
assumptions  which  are  valid  concerning  some  of  these  systems 
may  be  totally  unacceptable  when  referred  to  others.  For 
this  reason,  the  computer  operator  Is  given  a  considerable 
degree  of  freedom  to  choose,  by  the  way  In  iriilch  he  uses  the 
model,  the  assumptions  idilch  he  wishes  to  make  concerning 
the  Interrelation  of  various  system  parameters.  Nonetheless, 
there  are  some  assumptions  Inherent  In  the  development  from 
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which  there  Is  no  escape,  and  which  must  be  understood  If 
the  results  of  the  computer  program  are  to  be  Interpreted 
correctly.  Some  of  the  more  Important  of  these  assumptions 
are  listed  below. 

1 .  The  population  of  the  ground  state  Is  assumed  to  be 
eq.ual  at  all  times  to  the  doping  density;  that  Is, 
depopulation  of  the  ground  state  Is  negligible.  This 
assumption  Is  valid  as  long  as  the  doping  density  Is  much 
greater  than  the  Inversion  required  for  threshold. 

2.  The  entire  Laser  rod  Is  assumed  to  reach  threshold 
as  soon  as  the  average  population  Inversion  In  the  rod  Is 
sufflolent  to  sustain  the  dominant  mode.  The  validity  of 
this  assvimptlon  may  vary  from  one  material  or  geometry  to 
another,  and  should  be  remembered  \^en  restate  are 
Interpreted. 

3.  It  Is  assumed  that  the  pumping  pulse  can  be  repre¬ 
sented  In  time  by  the  sum  of  two  exponential  terms.  This 
assumption  was  made  after  studying  the  time  traces  of  xenon 
lamp  output  pulses,  and  It  was  concluded  that  a  time  function 
containing  two  exponential  terms  Is  adequate  for  the 
description  of  the  lamps  studied.  Its  acceptability  for 
other  pumping  sovircee  has  not  been  determined. 

4.  It  Is  assvuned  that  cai  coherent  output  Is  useful 
output,  and  no  attempt  Is  made  to  detezmlne  the  distribution 
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Of  energy  among  the  several  modes  i«hlch  may  be  exolted 


Approximations  In  the  Computer  Program 

Graphical  Integrations.  At  several  points  In  the 
computer  program  It  Is  necessary  to  evaluate  definite 
Integrals  by  graphical  (numerical)  means.  All  such 
Integrations  are  carried  out  by  using  the  trapezoidal  rule, 

I  ^f(x)dx  s  ^xf— f(xi+Ax)  +  f(x^+2nx) 

Jxi  2 

+  ...  ♦  f(x2-2nx)  +  f(x2-Ax)j  (3-1) 
with  one  exception,  a  case  where  tix^)  Is  Infinite.  In 
this  Instance,  the  Integral  Is  evaluated  as 

/>  x)dx  =^x[f(x^4Ax/2)  +  f(xi+3^/2) 

^  +  ...  +  f(x2-3£a/2)  +  f(x2-^x/2)]  (3-2) 

Although  there  are  certainly  more  accurate  graphical 
Integration  foimulas  than  these,  the  errors  arising  from 
Inaccurate  Integration  are  expected  to  be  minor  compared 
with  errors  In  the  experimental  determination  of  the  physical 
parameters  of  the  system  which  must  be  entered  as  Input 
data. 


Representation  of  Spectra.  Absorption  and  lamp  output 
spectra  are  represented  In  the  computer  program  by  the 
values  of  capture  cross-section  and  spectral  power  density 
at  Intervals  of  .01  microns  (100  i).  If  the  actual  spectra 
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being  represented  have  features  which  are  not  adequately 
described  by  such  a  representation,  some  adjustment  of  the 
Input  values  from  those  experimentally  determined  may  be 
necessary. 

Threshold  Time.  The  elapsed  time  between  the  beglzinlng 
of  the  pumping  pulse  and  the  achievement  of  threshold  enters 
into  the  calculation  of  radiation  efficiency.  In  the 
computer  program,  the  inversion  population  is  examined  for 
threshold  condition  at  intervals  of  20  microseconds,  and  the 
radiation  efficiency  is  calculated  using  threshold  time  to 
the  nearest  20  microseconds.  For  example,  if  threshold  is 
actually  achieved  in  164  microseconds,  the  computer 
calculates  the  efficiency  as  if  the  time  to  threshold  were 
180  microseconds,  since  that  is  the  first  time  checked  for 
which  the  inversion  population  is  sufficient  to  sustain  the 
dominant  mode. 

Undoubtedly,  a  considerable  improvement  in  accuracy 
could  be  effected  by  a  carefvil  revision  of  the  computer 
program  and  the  use  of  more  exacting  numerical  analysis.  It 
is  felt,  however,  that  the  program  as  it  stands  will 
represent  the  mathematical  model  adequately,  and  that  most 
of  the  errors  encountered  will  arise  from  the  Inadequacy  of 
the  model  itself.  Pexiiaps  the  least  certain  assumptions  of 
all  are  that  a  complex  mechanism  like  the  four-level  Laser 
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can  be  represented  by  anything  as  simple  as  the  model 
developed  In  this  paper,  and  that  the  physical  parameters  of 
the  system  can  be  measured  with  sufficient  accuracy  to  yield 
meaningful  results. 
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IV.  The  Efflolenoy  EQuatione 

In  Chapter  II,  the  overall  efficiency  of  a  four-level 
Laser  system  was  defined  as  the  product  of  six  sub¬ 
efficiencies.  The  object  of  this  chapter  is  to  develop 
mathematical  expressions  for  the  sub-efficiencies  in  terms 
of  measurable  physical  parameters  of  the  system. 

Independent  Variables  and  Constants 

Por  convenience,  the  system  parameters  will  be  grouped 
into  two  general  classes.  Those  parameters  which  are  more 
or  less  fixed  by  the  choice  of  Laser  material,  pumping 
method,  and  operating  environment  will  be  called  constants, 
and  those  idiich  are  more  or  less  under  the  control  of  the 
design  engineer  will  be  called  independent  variables.  From 
a  mathematical  standpoint,  there  is  very  little  difference 
between  these  two  classes  of  parameters;  the  main  difference 
is  reflected  in  the  computer  program  itself  by  the  fact  that 
once  the  constants  are  entered,  a  number  of  problems  can  be 
nm  with  different  values  of  the  Independent  variables 
without  reentering  the  constants.  This,  of  course,  does  not 
prevent  the  computer  operator  from  changing  constants 
whenever  he  feels  it  is  necessary. 

Independent  Variables.  The  following  system  parameters 
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have  been  selected  as  Independent  variables. 

1.  Laser  rod  length,  In  centimeters. 

2.  Laser  rod  radius,  R,  In  centimeters. 

3.  Active  Ion  density,  p.  In  Ions  per  cubic  centimeter. 

4.  End  mirror  reflectivities,  R^  and  R2»  In  watts  per 

watt. 

5.  End  mirror  absorptlvltles,  A^  and  A2>  In  watts  per 

watt. 

6.  Attenviatlon  constant  of  the  Laser  rod  material  at 
the  Laser  output  wavelength,  cxq.  In  nepers  per  centimeter. 

7.  Input  energy  to  the  pumping  sovirce,  E^ ,  In  joules. 

8.  Pvimplng  pulse  decay  time  constant,  In  seconds. 
Oonstants.  The  following  system  parameters  have  been 

selected  as  constants. 

1.  Lamp  efficiency,  as  defined  In  Chapter  II. 

2.  Geometrical  efficiency,  as  defined  In  Chapter 
II. 

3.  Lamp  spectrum,  p(^).  In  watts  per  centimeter. 

4.  Absorption  cross-section,  (7(A),  in  square 
centimeters. 

5.  Laser  output  wavelength,  A^,  In  centimeters. 

6.  Quauitum  efficiency,  &,  In  photons  per  photon. 

7.  Lower  euid  upper  limits  of  absorption  band,  emd 
A2,  In  centimeters. 
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8.  Energy  difference  between  ground  state  and  terminal 
state,  ^E,  in  Joules. 

9.  Lifetime  of  the  metastable  state,  "t*,  in  seconds. 

10.  Fluorescent  llnewldth  of  the  metastable -terminal 
transition,  hf,  in  cycles  per  second. 

11.  Ambient  operating  temperature,  T,  in  degrees  Kelvin. 

12.  If  the  lamp  spectrum  is  to  be  represented  by  black 
body  radiation,  the  temperature  of  the  black  body,  in 
degrees  Kelvin. 

It  may  seem  to  the  reader  that  listing  the  lamp 
efficiency  and  geometrical  efficiency  as  constants  is  an 
oversimplification.  Certainly  the  lamp  efficiency  will  vary 
with  the  energy  of  the  pulse,  and  the  geometrical  efficiency 
will  change  when  the  physical  dimensions  of  the  Laser  rod 
are  altered.  However,  while  these  two  quantities  are  not 
really  constant,  neither  is  their  relation  to  the  rest  of 
the  system  parameters  capable  of  being  expressed  in 
mathematical  terms.  If  the  actual  relationships  between 
these  quantities  and  the  other  parameters  are  known,  it  is 
possible  to  adjust  for  variations  in  lamp  efficiency  and 
geometrical  efficiency  by  applying  a  correction  factor  to 
the  value  of  lamp  input  energy  entered  into  the  program. 
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Capture  Efficiency 

Since  the  lamp  efficiency  and  the  geometrical 
efficiency  have  been  defined  as  constants,  the  first  step  In 
the  construction  of  the  model  Is  the  development  of  an 
expression  for  capture  efficiency.  This  will  first  be  done 
for  monochromatic  light,  calling  the  monochromatic  captvire 
efficiency  Wp,  and  then  extended  to  the  general  case 


Consider  a  Laser  rod  having  the  geometry  shown  In 
Figure  2.  Under  the  assxmptlon  (see  Chapter  III)  that  light 
entering  the  rod  Is  uniformly  spread  over  a  solid  angle  of 
arf  steradlans,  the  flux  arriving  at  point  p(r,0,0}  from  an 
element  of  surface  area  dA  located  at  point  q(R,0,z)  Is 

d'f'  =  e"^®  dA  watts/cm^  (4-1) 
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where  W  Is  the  total  Illuminating  power,  A  Is  the  total 
surface  area,  s  Is  the  distance  from  p  to  q,  and  Is  an 
attenuation  constant.  In  terms  of  the  Independent  variables 
and  constants  defined  earlier  In  this  chapter. 


A  =  2TTRi  cm^  (4-2) 

s  =  ( z2+r2+R2.2rR  cos  0)^  cm  (4-5) 

dA  =  R  dJ2(  dz  cm^  (4-4) 

o<.  =  ptf(A)  nepers/cm  (4-5) 


Substituting  (4-2),  (4-3) »  and  (4-4)  In  (4-1)  and 
integrating  ever  the  surface  gives  the  total  flux  at  the 
point  p  as 


''' 


■+Tr  |^^-oc(z2^.y2^R2.2rR  opg  0)k 

z2+r2+R2-2rR  cos  0 

2 

watts/cm 


dz  i0 


(4-6) 

Assuming  thatj^^^R,  or  that  both  ends  are  nearly  perfectly 
reflecting  surfaces,  and  taking  advantage  of  symmetry,  (4-6) 
can  be  written 


yO  -'O  z2+r2+R2-2rR  cos  0 

The  absorbed  power  density  at  any  point  In  the  rod  Is 
dW^  =  o(<f^  watts/cm^ 
and  the  total  power  absorbed  Is 


(4-7) 


(4-8) 


r  dr 


watts 


(4-9) 


« 
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The  monochromatic  captiire  efficiency,  by  definition,  Is 

Substituting  (4-7)  and  (4-9)  In  this  relation  gives 

Si' 


0  yo 


Equation  (4-10) 
variables.  Let 


^OO 


z^+r^+R^-2rR  cos  0)'^ 


dz  i0  dr  (4-10) 


yo  z2+r2+H2_2rR  cos  0 
can  be  simplified  considerably  by  change  of 


=  r/R 

(4-11) 

=  (r2+R2-2rR  cos  0)5 

(4-12) 

=  arctan  (z/a) 

(4-13) 

=  a/H  =  (1+R'2-2R'  cos  0)'^ 

(4-14) 

When  (4-11)  through  (4-14)  are  substituted  in  (4-10) 
and  the  results  simplified,  the  expression  for  '>1^  becomes 

^  (4.15) 


Equation  (4-15)  completely  specifies  Wp,  the  monochromatic 
capture  efficiency.  In  terms  of  one  parameter,  <xR,  Eq\iatlon 
(4-5)  shows,  however,  that  is  in  turn  a  fvinctlon  of  the 
wavelength  \  Hence  Wp  is  itself  a  function  of  A. 

If  the  spectral  power  density  of  the  pumping  source  Is 
p('A),  then  the  total  power  In  the  spectrum  between  the 
limits  A^  and  Ag  is 


W  = 


p(A)  dA 


watts 


(4-16) 
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The  power  absorbed  by  the  Laser  rod  is 


)A, 


V?r(^)  P(^)  watts 


(4-17) 


Further,  If  p(A)  Is  normalized  so  that  Its  average  value  over 
the  range  from  to  Ag  is  unity,  then  (4-16)  becomes 

W  =  A2-A^  watts  (4-18) 

and  the  capture  efficiency  is  given  by 

'Aa 


=  VW  =  ^  Wr(A)  p(A)  dA 


(4-19) 


Relaxation  Efficiency 

In  Chapter  II,  it  was  established  that  the  relaxaticn 
efficiency  is  equal  to  the  quantum  efficiency  £  multiplied 
by  the  ratio  of  the  energy  of  an  output  photon  to  the 
average  energy  of  the  pumping  photons.  Since  the  energy  of 
a  photon  is  inversely  proportional  to  its  wavelength,  the 
relaxation  efficiency  is  also  equal  to 

ri(j  =  (4-20) 

where  Ap  is  the  wavelength  of  a  photon  idiose  energy  is  the 
average  energy  of  the  pumping  photons. 

For  convenience,  let  the  unit  of  energy  in  the 
following  discussion  be  the  energy  associated  with  a  photon 
of  wavelength  one  centimeter  (1  om"^).  Then  the  photon 
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density  associated  with  a  spectral  power  density  p(A)  Is 

n(A)  = -^p(^)  photons/  cm'sec  (4-21) 

The  total  number  of  photons  absorbed  per  second  Is 
/'>2 

=  I  n(>v)  d)v  photons/sec  (4-22) 

and  the  total  power  absorbed  Is 

^a  ”  ^  cm"Vsec  (4-25) 

The  average  energy  per  photon  Is  Wg/N^,  and  the  wavelength 
of  the  corresponding  "average"  photon  Is  the  reciprocal  of 
the  energy, 

Ap  =  Na/*^a  “  (^  \  )V1  L  Ap(X)Wp(^)  dX  cm  (4-24) 

Since  the  capture  efficiency  appears  in  the  denominator  of 
the  expression  for  Xp,  It  will  also  appear  In  the 
denominator  of  the  expression  for  relaxation  efficiency. 
Further,  since  the  model  makes  no  direct  use  of  either  T|q  or 
n,Q  Individually,  considerable  time  and  effort  can  be  saved 
by  computing  the  product  directly.  Combining  equations 

(4-20)  and  (4-24), 

/a"' 

Radiation  Efficiency 

The  radiation  efficiency  was  defined  in  Chapter  II  as 
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the  ratio  of  the  energy  which  enters  the  coherent  field  to 
the  total  energy  involved  in  the  transition  from  the 
metastable  state  to  the  terminal  state.  Since  all  photons 
emitted  by  ions  making  this  transition  are  of  the  same 
energy,  the  radiation  efficiency  is  also  equal  to  the 
ratio  of  the  number  of  photons  contributing  to  the  coherent 
field  to  the  total  niuaber  emitted  at  the  Laser  output 
wavelength. 

When  a  pulse  of  energy  from  the  pumping  source  is 
applied  to  the  four-level  Laser,  the  population  of  the 
metastable  state  increases  with  time  until,  if  the  pumping 
energy  is  sufficiently  high,  a  point  is  reached  >Aiere  the 
gain  of  the  system  for  the  coherent  field  exceeds  unity. 
Luring  the  buildup  of  population,  some  ions  are  also 
decaying  spontaneously  to  the  terminal  state,  and  the 
energy  associated  with  the  transition  of  these  ions  is 
lost.  After  the  threshold  has  been  reached,  however,  the 
coherent  field  builds  up  vezy  rapidly,  and  nearly  all  of 
the  ions  idxich  make  the  transition  thereafter  contribute 
their  energy  to  the  coherent  field. 

To  determine  the  radiation  efficiency,  then,  it  is 
necessary  to  (1)  obtain  an  expression  for  the  population 
of  the  metastable  state  as  a  function  of  time,  (2)  determine 
the  population  required  for  threshold,  (?)  find  the  time  for 
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which  the  two  are  eq.ual,  and  (4)  calculate  the  amount  of 
energy  lost  prior  to  threshold  and  thence  obtain  the 
radiation  efficiency. 

Metastable  State  Population.  Assume  (see  Chapter  III) 
that  the  pumping  pulse  can  be  represented  by  the  siun  of  two 
exponential  terms, 

W(t)  =  -  e“'‘^/'^2)  watts  (4-26) 

The  constant  can  be  foimd  by  imposing  the  condition 


rw(t)  dt  =  Joules 


(4-27) 


When  this  Is  done,  the  value  of  Is  obtained. 


(4-28) 


The  equation  for  the  population  of  the  metastable  state  Is 


||  =  KW(t)  -  n/r  lons/sec  (4-29) 


where 


K  s  lons/watt ‘sec 

he 


(4-30) 


If  (4-26),  (4-28),  and  (4-30)  are  substituted  In  (4-29),  the 
resulting  equation  can  be  solved  for  n,  giving 

E  Ao  T2e”'^A'2 

“  “  110^1-2^)  t  (T,.r)  *  (Tg-T) 
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Since  this  equation  represents  the  sum  total  of  the 
population  of  the  metastable  state,  the  average  density  of 
Ions  In  the  metastable  state  Is  n/V,  where  V  Is  the  voltime 
of  the  rod  In  cubic  centimeters. 

Threshold  Condition.  All  of  the  preceding  development 
applies  equally  well  idiether  the  Laser  system  Is  used  as 
an  oscillator,  providing  Its  own  excitation,  or  as  an 
amplifier,  to  build  up  an  existing  coherent  field.  The 
threshold  conditions,  however,  are  different  for  these  two 
modes  of  operation,  and  since  It  Is  desired  that  the  model 
represent  both  the  Laser  oscillator  and  the  Laser  amplifier, 
both  threshold  conditions  will  be  developed. 

When  an  electromagnetic  wave  of  the  proper  frequency 
Is  passed  through  the  Laser  material.  If  the  population  of 
the  metastable  state  exceeds  the  population  of  the  terminal 
state,  the  wave  will  be  amplified  by  the  stimulated  emission 
of  radiation  according  to  the  equation 

W(z)  =  (4-32) 

where  z  Is  the  direction  of  propagation,  Is  the 
attenuation  constant  of  the  material,  and 

y  =  nepers/cm  (4-33) 

In  (4-33) >  Is  the  population  difference  between  the 
metastable  state  and  the  terminal  state  In  Ions  per  cubic 
centimeter,  and  ^^2  ^s  the  capture  cross-section  between 
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the  two  states,  given  by  (Ref  4:A-2) 

Applying  equation  (4-32)  to  the  Laser  oscillator  of  Figure 
3,  the  following  relations  are  obtained: 


W4  =  (4-37) 

If  the  gain  of  the  system  is  unity,  this  also  means  that 

=  R2^^4  =  (4-38) 


The  solution  of  (4-38)  for  Y  gives 
y  20(qS.  -  ln(RiR2) 

2  =  - 2T - 


(4-39) 


whence,  by  equation  (4-33), 


= 


20fai?  -  ln(RiR2) 


2^2^ 


(4-40) 


The  total  ion  density  In  the  metastable  state  required  for 
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threshold  Is  obtained  by  adding  to  An  the  thermal  population 
of  the  terminal  state,  thus 

at  =  p.-“Al  *  £2aLl^^lM  lons/omJ  (4-41) 

The  threshold  condition  can  be  found  from  equations  (4-31 ) 
and  (4-41):  that  Is,  threshold  Is  achieved  uhen 

^  =  at  (^--^2) 

Threshold  for  a  Laser  amplifier  will  be  defined  as  the 
condition  such  that  the  net  gain  of  the  amplifier  Is  unity. 
Neglecting  reflections  at  the  end  surfaces,  and  applying 
equation  (4-32), 

=  ^<0  t^-^3) 

n^  =  +  Ob/<^2  lons/cm^  (4-44) 

Again,  equation  (4-42)  defines  the  required  condition  for 
threshold. 

Threshold  Time.  The  solution  of  equation  (4-42) 
explicitly  for  time  Is  a  problem  of  considerable  magnitude. 
Although  this  solution  could  be  obtained  by  Iteration  to  any 
desired  accuracy,  some  compromise  Is  necessary  between 
obtaining  an  acciirate  value  for  threshold  time  and  keeping 
the  time  required  for  calculation  down  to  a  reasonable 
level.  As  discussed  In  Chapter  III,  the  actual  procedure 
selected  Is  to  Insert  values  of  time  at  20  mloroseoond 
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Intervals  and  check  to  see  If  threshold  has  been  achieved. 
This  process  will  be  discussed  more  fully  In  Chapter  V. 


1 


Computation  of  Radiation  Efficiency.  If  the  time  to 
threshold  Is  found  to  be  seconds  for  a  particular  Laser 
system,  the  number  of  Ions  which  will  have  decayed 
spontaneously  prior  to  Is  given  by 


n/T  dt 


Ions 


(4-45) 


Substituting  the  value  of  n  from  equation  (4-31)  and 
performing  the  Integration  gives 


The  total  nimber  of  Ions  making  the  transition  Is 


,  EWjsiyia??  ion. 
he 


and  the  radiation  efficiency  Is  given  by 

=  (Ni  -  Nl)/Ni  (4-48) 

Substitution  of  (4-46)  emd  (4-4?)  in  (4-48)  gives  finally 


tir  .  r2r] 


(1 


(4-49) 
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Equation  (4-49)  represents  the  radiation  efflolency  of 
a  Laser  oscillator;  however,  the  radiation  efficiency  of  a 
Laser  amplifier  will  depend  upon  how  much  excitation  Is 
supplied  at  the  Input.  If  the  exciting  signal  Is  very 
large,  the  efficiency  will  be  fairly  well  represented  by 
equation  (4-49),  and  the  gain  will  be  very  low.  On  the 
other  hand,  If  the  exciting  signal  Is  very  small,  the 
efficiency  will  be  very  low,  and  the  gain  will  be  higher. 

An  Investigation  of  the  relation  between  gain  and  efflolency 
In  the  Laser  amplifier  Is  beyond  the  scope  of  this  study. 
However,  so  that  the  model  will  be  of  some  use  to  the 
designer  of  Laser  amplifiers,  two  points  on  the  gain- 
efficiency  curve  will  be  calculated,  xiamely,  the  efficiency 
for  unity  gain  and  the  gain  for  zero  efficiency.  These  are 
labelled  the  maximum  efflolency  and  maximum  gain 
respectively,  and  It  must  be  emphasized  that  the  two  cannot 
exist  simultaneously. 

The  gain  for  very  small  signals,  from  equation  (4-32), 
is 

®max  =  e^‘^32^%ax-‘^o)^  (4.50) 

where  An^^x  Is  the  maximum  inversion  achieved  during  the 
pumping  period. 
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Optical  Efficiency 

Vfith  five  cf  the  six  sub-efficiencies  defined,  either 
as  ccnstants  or  in  terms  of  the  system  parameters,  the  only 
step  remaining  in  the  mathematical  development  of  the 
efficiency  equations  is  the  derivation  of  an  expression  for 
the  optical  efficiency. 

The  optical  losses  in  the  system  include  both 
attenuation  by  the  Laser  rod  Itself  and  reflector  losses  at 
the  ends  of  the  rod  (if  the  Laser  is  used  as  an  oscillator). 
Referring  again  to  Figure  3*  page  25,  and  for  convenience 
making  the  substitution 

g  =  y  -  %  (^-51 ) 

the  attenuation  losses  in  the  leftward  travelling  wave  aoid 
the  rightward  travelling  wave  are  respectively 


*^AR  =  Jq 


.gx 


dx 


Wje®*  dx 


(4-52) 


(4-53) 


When  equaticns  (4-36)  and  (4-51)  are  combined  with  (4-52) 
and  (4-53)  and  the  resulting  equations  integrated,  the  total 
loss  due  to  attenuation  is  found  to  be 


^LA  =  ^AL*^JLR  =  ^  ^  O  ( 4-54) 

The  reflector  losses  at  the  left  end  and  right  end 
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Wfil,  =  A,  *2  (A-55) 

Waa  =  *2''* 

Substitution  of  (4-35),  (4-36),  and  (4-51)  Into  the  above 
gives  the  total  reflector  loss  as 

WiR  =  ifiiL  ♦  Was  =  +  Aa/Ba)  (4-57) 

The  total  loss  Is 

Wi  =  Wla  ♦  Bib 

and  the  optloal  efflolencir  Is 

\  =  tfo/(Wo  +  Wl)  (4-59) 

where  Wq  Is  the  output  power,  given  by 

Wq  =  W4(1-R2-A2)  =  (I-R2-A2)  (4-60) 


When  (4-54)  through  (4-60)  are  combined  and  simplified,  the 
resulting  expression  for  optical  efficiency  becomes 


n-  _ _ y  (4-61) 

"  g(1-B2)  V  <3oB2(1eBieEe)(eE‘-1)  e  AiBaeES 


The  value  of  g  can  be  obtained  from  equations  (4-59)  and 


(4-51)  as 


_  -ln(RtR2) 

52 


(4-62) 


Since  the  maximw  efficiency  of  the  Laser  amplifier  is 
calculated  with  xinity  gain  assumed,  the  power  W  is 
independent  of  distsuice  along  the  rod,  and  the  total  loss 
(again  neglect ing  refleotions  at  the  ends)  is 


30 


GE/EE/62-18 

Wl  =  (4-63) 

The  optical  efficiency  of  the  Laser  amplifier  Is  then 

^0  =  V(Wo  +«foi^Wo)  =  1/(1  +^0^)  (^-64) 

The  so-called  "walkoff"  loss  has  been  neglected  In  the 
derivation  of  the  optical  efficiency.  Fox  and  LI  (Ref 
3:484)  have  shown  that  the  wallcoff  lose  for  the  dominant 
modes  In  a  Fabry-Perot  resonator  Is  negligible  compared 
with  the  attenuation  and  reflector  losses.  In  the  Laser 
amplifier,  the  severity  of  wallcoff  loss  will  be  determined 
by  the  beam  divergence  of  the  exciting  signal.  If  the  Input 
to  the  Laser  amplifier  Is  the  output  of  a  Laser  oscillator, 
the  walkoff  loss  will  be  very  slight. 

Summary 

In  this  chapter,  a  mathematical  model  for  the 
determination  of  Laser  efficiency  has  been  developed.  The 
model  consists  of  expressions  for  four  of  the  six  sub- 
efflclencles  which  were  defined  In  Chapter  II,  In  terms  of 
the  measurable  physical  parameters  of  the  Laser  system. 

These  sub-efficlencles,  when  multiplied  together  and  by  the 
lamp  efficiency  and  geometrical  efficiency,  which  have  been 
defined  as  constants,  give  the  overall  system  efficiency  of 
the  four-level  Laser  system. 
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V.  The  Oompttter  Progran 

The  mathematical  development  of  the  Laser  system  model 
having  been  completed,  the  next  step  Is  the  translation  of 
the  model  Into  a  program  for  the  IBM  1620  computer. 

Ob.lectlves 

In  the  formulation  of  the  computer  program,  It  has  been 
attempted  to  achieve  certain  objectives  and  conform  to 
oertaln  criteria  In  order  to  make  the  end  product  as  useful 
as  possible.  Specifically,  attainment  of  the  following 
objectives  was  attempted: 

1 .  The  final  program  should  be  capable  of  being  run 
on  any  IBM  1620  computer  with  card  Input/output,  with  no 
optional  features  required. 

2.  The  program  should  represent  as  closely  as  possible 
the  mathematical  model  developed  In  Chapter  IV. 

3.  Execution  of  the  program  should  be  as  convenient 
as  possible.  Ilherever  practicable,  the  program  should  be 
made  self-explanatory  by  the  use  of  typed  out  messages. 

4.  The  program  should  produce  results  fast  enough  to 
be  useful  as  a  design  aid;  that  Is,  during  one  reasonable 
period  of  operation.  It  should  be  possible  to  calculate 

the  efficiency  of  several  combinations  of  system  parameters. 
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Language 


The  program  has  been  written  in  a  modification  of  the 
IBM  FORTHM  language  devised  by  1st  Lt.  Pratt  of  the  Air 
Force  Institute  of  Technology  Mathematics  Department. 

There  are  a  few  differences  between  this  language  and 
standard  FORTRAN,  the  more  Important  from  the  programmer's 
standpoint  being  the  following: 

1 .  Format  specification  Is  not  required  on  Input 
statements,  unless  It  Is  desired  to  prevent  the  reading  of 
some  of  the  data  on  a  card. 

2.  The  automatic  typewriter  carriage  return  at  the 
beginning  of  each  PRINT  statement  has  been  deleted,  allowing 
a  full  87  character  line  to  be  typed  by  two  successive 
PRINT  statements. 

3*  Trace  Instructions  are  automatically  compiled,  and 
the  resulting  object  program  Is  at  least  as  short  as  the 
same  program  compiled  In  standard  FORTRAN  without  trace 
Instructions. 

4.  The  subroutines  have  been  shortened  so  that  the 
object  program  begins  at  storage  location  OSOOO. 

These  differences  have  been  noted  here  to  Indicate  that 
some  changes  In  the  source  program  statements  may  be 
required  If  It  Is  desired  to  compile  the  program  developed 
In  this  chapter  In  standard  FORTRAN. 
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lo  meet  the  first  objective  set  forth  above,  the 
program  must  be  capable  of  being  run  on  a  computer  having 
a  20,000  digit  memory.  lo  meet  this  objective.  It  was 
necessary  to  spilt  the  program  Into  three  smaller  programs 
which  are  run  consecutively.  The  divisions  have  been 
arranged  so  that  part  1  of  the  program  Is  completely 
Independent  of  the  system  parameters,  and  need  only  be  run 
once.  The  output  of  part  2  Is  dependent  only  on  certain  of 
the  system  parameters  ^Ich  have  been  defined  as  constants, 
and  this  part  need  only  be  rerun  when  the  constants  Involved 
are  changed,  fart  3  Is  that  portion  of  the  program  which 
accepts  the  values  of  the  Independent  variables  and 
calculates  the  final  results. 

The  intermediate  data  represented  by  the  outputs  of 
part  1  and  part  2  of  the  program  are  stored  on  cards  and 
read  Into  the  computer  at  the  beginning  of  execution  of  the 
following  part.  Thus  If  It  Is  desired  to  Investigate  the 
effect  of  changing  certain  variables  In  a  system,  and  part 
2  of  the  program  has,  at  some  previous  time,  been  run  using 
constants  appropriate  to  the  Laser  material  Involved,  It  Is 
only  necessary  to  read  the  appropriate  card  deck  Into  part  3 
and  enter  the  desired  values  of  the  variables. 
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Block  Diagram  of  Bart  1  of  Computer  Brogram. 
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Part  i 

Part  1  of  the  computer  program  couslets  of  the 
necessary  Instructions  to  effect  the  solution  of  equation 
(4-15), 

^  r  4K  4E'  (4-15) 

A  block  diagram  of  the  program  Is  shovn  In  Figure  4.  The 
program  first  constructs  a  table  of  the  values  of  the  Inner 
Integral  for  values  of  the  parameter  o(Ra' .  Integration  Is 
done  by  the  trapezoidal  rule,  equation  (3-1),  and  the 
Interval  AJT  used  Is  '^100.  The  Integral  Is  computed  for 
values  of  oiRa'  from  zero  to  ten  at  Intervals  of  0.05,  and 
the  resulting  values  are  stored  as  a  subscripted  variable. 
The  next  Integration  Is  carried  out  with  o(E  and  R'  as 
parameters.  Integration  Is  accomplished  uslzxg  equation 
(3-2)  (rectangular  elements  of  area),  and  ^  Is  ^/lOO.  R' 

Is  allowed  to  vary  from  zero  to  one  In  steps  of  0.05,  and  «R 
takes  on  values  from  0.01  through  10.24,  each  of  which 
(except  the  first)  Is  V?  times  the  preceding  value.  The 
last  Integration  uses  the  trapezoidal  rule,  being  0.05. 
The  output  consists  of  a  set  of  values  of  one  for  each 
value  of  cxR,  punched  one  to  a  oard  (21  cards).  The  actual 
source  statements  used  in  the  compilation  of  this  and  the 
other  two  parts  of  the  program  are  listed  In  Appendix  A. 
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Part  1  of  the  program  req.uires  13  hours  for  execution, 

and  since  It  has  already  been  run,  there  Is  no  need  to 

run  It  again  unless  It  Is  desired  to  alter  the  logic  In 

some  vay.  The  required  cards  for  Input  to  part  2  may  be 

punched  manually,  and  If  this  Is  done,  the  cards  should  be 

punched  in  colxmms  2  through  14  as  follows: 

.1 93971 79E-01 

. 272096 17E-01 

.38035933E-01 

.529061803-01 

.731800463-01 

.10051629E+00 

. 1 3685365E+OO 

.184181663+00 

.24438795E+00 

.318499943+00 

.4059441 5E+00 

.50361047E+00 

.605315163+00 

.70225637E+00 

.785719703+00 

.849706043+00 

.893819303+00 

.921 8641 5E+00 

.939198063+00 

.947537623+00 

.948452263+00 

The  monochromatic  capture  efficiency  ranges  from 
about  2  percent  for  aR  =  0.01  to  about  95  percent  for 
o(R  =  10.24.  It  Is  believed  that  this  range  suffices  to 
represent  most  commonly  used  four-level  Laser  systems. 

If  It  Is  desired  to  run  this  part  of  the  program  In 
the  future,  all  that  Is  required  on  the  part  of  the  operator 
Is  to  load  the  program  and  press  PUIICH  START  and  START. 
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Fig.  5 

Slook  Slagran  of  Part  2  of  Computer  Program. 


GE/BE/62-18 


During  execution  of  the  prograffl*  several  error  messages 
will  be  typed  out,  such  as  ERROR  E4  and  ERROR  F5*  These 
messages  arise  from  computations  of  the  form 


y  = 


,-(1/0) 


(5-1) 


and  should  be  Ignored,  as  the  correct  value  of  y,  namely, 
zero,  results  from  this  computation. 


Part  2 

Part  2  of  the  computer  program  accepts  the  output 
cards  from  part  1  and  some  of  the  system  parameters  ^ioh 
were  defined  as  constants  In  Chapter  IV,  axxd  computes  the 
solution  of  equation  (4-25) » 

°  (*-25) 

1  block  diagram  of  the  program  Is  shown  In  Figure  5*  After 
accepting  via  the  typewriter  the  values  of  &,  ^2* 

the  program  determines,  from  the  setting  of  a  console 
switch,  whether  the  operator  wishes  to  enter  the  lamp 
spectrum  explicitly  or  allow  It  to  be  represented  by  black 
body  radiation.  In  the  former  case,  the  program  reads  the 
values  of  (T  and  p  for  the  range  from  \  to  f^om  cards 
irtilch  have  been  manually  punched  with  the  appropriate  values. 
It  then  nozmallzes  p(A)  so  that  Its  average  value  is  unity. 

In  the  latter  case,  ozily  is  read  from  ceirds,  and  the 
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computer  generates  the  black  body  radiation  function 
“  e{hcAT^)-1 

at  a  temperature  selected  by  the  operator,  where  K  Is 
adjusted  so  that  the  average  value  of  p(X)  is  unity. 

Kext,  the  program  reads  the  output  cards  from  part  1 
and  Interpolates  the  table  of  to  obtain  the  correct  value 
for  each  wavelength,  with  as  a  parameter.  pR  Is  allowed 
to  vary  from  0.001x10^®  to  1,024x10^®  cm“^,  with  each  value 
(except  the  first)  being  V?  times  the  preceding  value.  The 
range  of  <r  which  the  program  will  handle  without  distortion 
varies  with />R.  Values  which  result  In  o^R  less  than  0.01 
are  assigned  =  0,  and  those  which  result  In  values  of  o(B. 
greater  than  10.24  are  assigned  the  highest  value  In  the 
table,  Wy  =  0.94845226.  The  particular  range  selected  was 
designed  to  fit  the  trlvalent  neodymium  Ion,  iihlch  has  a 
peak  absorption  cross-section  of  about  8x10"^®  cm^  (Ref 
1:196). 

The  Integration  Is  performed  by  the  trapezoidal  rule 
with  .01  micron.  The  output  consists  of  a  value  of 

for  each  value  of  pR,  punched  one  to  a  card  (21  cards). 
Operation.  The  operation  of  part  2  of  the  program  Is 
as  follows:  jLfter  the  program  Is  loaded  and  the  START  button 
pressed,  the  following  message  Is  typed: 
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SWITCH  1  011  TO  HEPSESENT  LAI'i?  SPECTRUM  BY  BLACK  BODY 
EADUTIOII,  PUSH  START. 

and  the  computer  halts.  After  this  Instruction  has  been 
compiled  with,  the  following  message  Is  typed: 

ENTER  QUANTUM  EPPICIENCY. 

The  quantum  efficiency  of  the  Laser  material  Is  entered  as  a 
decimal  fraction.  The  following  message  Is  then  typed: 

ENTER  LH'IITS  OF  ABSORPTION  SPECTRUM  IN  MICRONS. 

The  limits  and  ^2  are  entered;  the  order  of  their  entry 
Is  Immaterial.  The  following  constraints  are  placed  on  the 
values  of  these  limits: 

1 .  ^2  ~  must  be  an  Integral  multiple  of  .01 
microns. 

2.  ^2  **  must  not  exceed  1  micron. 

After  these  values  have  been  entered,  the  following  message 
Is  typed: 

ENTER  LASER  OUTPUT  WAVELENGTH  IN  MICRONS. 

Xq  Is  then  entered. 

At  this  point.  If  console  switch  1  Is  on,  the  following 
message  Is  typed: 

EilTER  BLACK  BODY  lEl-IPSRAIURE  IN  DEGREES  K. 

Tjj  Is  entered,  and  the  black  body  radiation  function  Is 
generated.  The  program  will  then  read  from  cards  the 
values  of  a(X) .  The  Input  cards  must  contain  ONLY  the  values 
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— pn  2 

of  ^  in  units  of  10  cm  ,  punched  one  to  a  card,  and 
arranged  In  order  of  Increasing  wavelength  from  X|  to 
inclusive.  Each  card  read  will  be  associated  with  a  wave¬ 
length  0.01  microns  longer  than  the  previous  card,  and  the 

program  will  read  exactly  +  1  cards. 

.01 

If  console  switch  1  is  off,  the  program  will  read  from 
cards  the  values  of  and  p(>v).  The  cards  must  be 

punched  with  two  numbers  each,  the  first  number  being  6^ 
in  10  ^  cm  and  the  second  being  the  spectral  power  density 
p  in  any  units  desired  (the  computer  will  normalize  the 
spectrum) .  The  vravelength  value  associated  with  each  card 
and  the  order  of  the  cards  is  as  described  in  the  preceding 
paragraph. 

At  this  point,  the  output  cards  from  part  1  will  be 
read,  and  the  remainder  of  the  program  will  proceed  without 
further  intervention  by  the  operator.  The  total  running 
time  for  this  part  of  the  program  (using  an  absorption  band 
0.6  microns  wide)  is  about  25  minutes. 

x'art 

jrart  3  of  the  program  accepts  the  output  data  from 
part  2  and  the  remaining  constants  and  Independent  variables, 
and  computes  either  the  efficiency  of  a  Laser  oscillator  or 
the  maximum  efficiency  and  maximum  gain  of  a  Laser  amplifier. 
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depending  on  the  setting  of  a  console  switch.  A  block 
diagram  of  part  3  is  shown  In  Figure  6. 

At  the  beginning  of  the  program,  the  output  cards  from 
part  2  are  read,  the  constants  >q, 

T  are  accepted,  and  some  Intermediate  constants  are 
computed  to  save  time  later  In  the  program.  The  variables 
<%and  p  are  accepted,  and  the  thermal  population  of  the 


terminal  state. 


Pt  = 


(5-3) 


Is  computed.  The  variables  R  and  2  are  accepted,  and  a 
check  Is  made  to  see  if  .001  ^pEbclO"^®^  1 .024.  If  not, 
new  values  of  p,  R,  and  £  are  requested. 

At  this  point.  If  console  switch  1  Is  on,  the 
population  required  for  threshold  Is  computed  from  equations 
(4-44)  and  (5-3)  as 

nt  =  Pt  +  V(?32  ^5-4) 

If  the  switch  Is  off,  R2,  A^ ,  and  A2  are  accepted.  R]  is 
assumed  to  be  1  -  A^  (no  transmission  through  one  end) ,  and 
the  threshold  population  Is  calculated  using  equations 


(4-41)  and  (5-3), 


at  =  Pt  + 


2oro£  -  ln(R,R2) 


(5-5) 


The  values  of  E^  and  are  then  accepted.  The  value 
of  'i'2  calculated  as  /1 0.35*  This  ratio  was  arrived  at 
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on  the  basis  of  oscilloscope  traces  of  the  output  of  a 
General  Electric  FT-506  helical  xenon  flash  lamp,  and  Is 
the  ratio  of  time  constants  which  produces  a  power  maximum 
at  time  t  =  .  3ome  of  the  oscilloscope  photographs 

used  as  the  basis  for  the  selection  of  the  time  constant 
ratio  are  shown  in  Figure  7.  The  traces  shown  are  all  of 
50c  ;]oule  pulses,  the  Input 
voltages  and  capacitances 
being,  from  top  to  bottom, 

(1  )  1225  V,  672  /If. 

(2)  1000  V,  1008 /If. 

(3)  567  V,  1344  /If. 

(4)  735  V,  1680  Mf. 
x'he  time  scale  is  1  cm=.5  ms. 

If  for  any  reason  It  Is 
desired  to  change  the  ratio, 

It  is  only  necessary  to 
reflect  the  cheuige  In 
statement  number  120  in  the 
source  program  (see  Appendix 
A) . 

it  Is  interesting  to 

note  that  for  the  lamp  checked,  In  microseconds  is 
approximately  equal  to  the  capacitance  in  microfarads. 


Fig.  7 


Typical  Output  Waveforms 
of  the  G.E.  FT-506  Lamp. 
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ilavlnf?  obtained  the  values  of  all  of  the  system 
parameters,  the  program  Interpolates  the  table  of  to 

obtain  the  correct  value  for  the  radius  and  Ion  density 
entered,  and  proceeds  to  calculate  the  metastable  state 
population  at  Intervals  of  20  microseconds,  comparing  each 
result  with  the  value  of  If  threshold  Is  not  reached 

In  640  microseconds,  a  message  Is  typed  Indicating  that 
more  energy  Is  required.  The  time  for  which  the  value  of 
the  metastable  state  population  (n/V)  first  exceeds  or 
equals  n.^^  Is  considered  to  be  ,  the  threshold  time,  and 
either  the  efficiency,  equations  (4-49),  (4-6l),  and  (2-1), 
or  the  maximum  efficiency  and  maximum  gain,  equations  (4-49), 
(4-64),  (2-1),  and  (4-50),  are  computed,  depending  on  the 
setting  of  console  switch  1 ,  In  the  computation  of  maximum 
gain,  the  metastable  population  Is  computed  at  20  microsecond 
Intervals  until  a  value  Is  found  which  Is  less  than  the 
previous  value.  The  higher  value  Is  then  used  to  compute  the 
maximum  gain.  If  the  population  is  still  Increasing  at 
t  =  640  microseconds,  the  value  of  n  corresponding  to 
t  =  640  microseconds  is  used  in  the  computation. 

After  typing  the  results,  the  computer  will  halt. 
Pressing  the  STAHT  button  returns  the  program  to  the  point 
where  oCq  was  entered,  and  the  cycle  repeats. 

Operation.  The  operation  of  part  5  of  the  program  Is 
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as  follows:  After  the  program  has  been  loaded  and  the  START 

button  pressed,  the  output  cards  from  part  2  are  read  and 

the  following  message  is  typed: 

ENTER  ENERGY  DIPPERENOE  BETWEEN  GROUND  STATE  AIJD 
TERMINAL  STATE  IN  RECIPROCAL  CM. 

The  value  of  hB  Is  entered,  and  the  following  message  Is 

typed: 

ENTER  LASER  OUTPUT  WAVELENGTH  IN  MICRONS. 

Is  entered,  and  the  following  message  Is  typed: 

ENTER  LIPETIME  OP  METASTABLE  STATE  IN  MICROSECONDS. 

Is  entered,  and  the  following  message  Is  typed: 

ENTER  LAMP  EPPICIENCY,  GEOMETRICAL  EPPICIENCY. 
and  YIq  are  entered.  It  should  be  remembered  that  In  the 
determination  of  only  the  portion  of  the  spectnun 
between  \  and  ^2  Is  considered  as  useful  output.  After 
these  values  have  been  entered,  the  following  message  Is 
typed: 

ENTER  PLUORESOENT  LINE  WIDTH  IN  CYCLES  PER  SBCOIUD. 

Af  Is  entered,  and  the  following  message  Is  typed: 

ENTER  OPERATING  TE2^IPERATURE  IN  DEGREES  K. 

T  Is  entered,  and  the  following  message  Is  typed: 

ENTER  ATTBIJUATION  CONSTANT  OP  HOST  MATERIAL. 

(Xq,  In  nepers/cm.  Is  entered,  and  the  following  message  Is 
printed : 
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ENTER  DOPING  DENSITY  IN  IONS  PER  CO. 
jO  is  entered,  and  the  following  message  is  typed: 

ENTER  ROD  RADIUS  AND  LENGTH  IN  OM. 

R  and  H  are  entered,  and  a  oheok  is  made  to  see  if  pR  is 
within  the  acceptable  range.  If  not,  the  following  message 
is  typed: 

RHO  R  OUT  OP  RAi^GS 

ENTER  DOPING  DENSITY  IN  IONS  PER  00. 

As  soon  as  an  acceptable  combination  of  p  and  R  has  been 

obtained,  the  following  message  is  typed: 

SWITCH  1  ON  POR  A>IPLIPIER,  OPP  POR  OSCILLATOR, 

PUSH  START. 

If  console  switch  1  is  set  to  the  ON  position,  the  remainder 
of  the  program  will  calculate  the  mazimiun  efficiency  and 
maximim  gain  of  a  Laser  amplifier.  If  the  switch  is  set  to 
OPP,  the  program  will  calculate  the  efficiency  of  a  Laser 
oscillator. 

After  the  switch  has  been  set  and  the  START  button 
depressed,  if  switch  1  is  off,  the  following  message  is 
typed: 

ENTER  RBPLEOTIVITY  OP  OUTPUT  END. 

R2  is  entered,  and  the  following  message  is  typed: 

ENTER  ABSORPTIVITY  OP  CLOSED  END,  OUTPUT  END. 

A)  and  entered.  If  switch  1  is  on,  the  reflectivities 

and  absorptivities  are  not  required  by  the  program  (assumed 
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zero),  and  the  preceding  two  messages  are  omitted.  In  any 

event,  the  following  message  Is  then  typed: 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TH-IE  CONSTANT 

IN  MICROSECOI'IDS. 

E^  and  are  entered,  the  computer  proceeds  to  compute  the 

metastable  population,  and  the  results  are  compared  with  the 

threshold  condition.  If  threshold  Is  not  reached,  the 

following  message  Is  typed: 

THRESHOLD  NOT  REACHED,  INCREASE  LAI'^IP  ENERGY. 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT 

IN  MICROSECONDS. 

When  a  value  of  lamp  energy  has  been  obtained  which  produces 
the  threshold  condition,  the  following  message  Is  typed: 

I  DIE  TO  THRESHOLD  =  XXX  MICROSECONDS 
The  efficiency  (switch  1  off)  or  maximum  efficiency  and  gain 
(switch  1  on)  are  computed,  and  the  results  are  printed  In 
the  appropriate  one  of  the  following  two  formats: 

EFFICIENCY  s  X.XXXXXXXX  PERCENT. 

MAX.  EFFICIENCY  =  X.XXXXXXXX  PERCENT,  MAX.  GAIN  = 

xxx.xxxxxxx 

After  printing  the  results,  the  computer  halts.  A  new  set 
of  Independent  variables  may  be  entered  by  pressing  the 
START  button.  The  running  time  for  part  3  Is  about  eight 
minutes  for  the  first  set  of  parameters,  and  two  to  three 
minutes  for  each  subsequent  set  of  independent  variables. 
Including  time  required  to  manually  enter  the  required  data. 
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Summary 

In  this  chapter,  a  computer  program  has  been  described 
which  calculates  the  efficiency  of  a  Laser  oscillator  or 
the  maximum  efficiency  and  maximum  gain  of  a  Laser  amplifier 
from  the  values  of  the  system  parameters.  The  program  can 
be  run  on  any  IBli  1620  computer  which  has  card  input/output, 
and  follows  olosely  the  mathematical  model  developed  in 
Chapter  IV.  It  is  convenient  to  use,  and  has  been  made 
almost  completely  self-explanatory  by  the  use  of  typed  out 
messages.  After  about  30  minutes  of  precomputation,  the 
program  will  process  a  set  of  independent  variables  and 
print  the  results  in  two  to  three  minutes. 
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VI.  A  Sample  Problem 

This  chapter  Is  an  Illustration  of  how  the  computer 
model  developed  In  the  preceding  chapters  might  be  applied 
to  a  design  problem.  It  Is  not  Intended  to  be  an  actual 
analysis  of  any  specific  Laser  system,  existing  or  proposed, 
although  an  effort  has  been  made  to  make  the  values  of  the 
constants  used  In  the  sample  problem  conform  as  closely  as 
possible  to  the  known  characteristics  of  the  materials 
Involved.  In  some  cases,  these  constants  were  actually 
measured  by  the  author.  In  other  cases  they  were  extracted 
from  the  literature,  and  In  still  others,'  they  represent 
pure  guesswork. 

Nonetheless,  It  Is  believed  that  all  of  the  system 
parameters  used  In  the  model  are  physically  measurable,  and 
It  would  be  Interesting  to  actually  measure  these  quantities 
for  a  real  Laser  system  and  compare  the  measured  efficiency 
under  various  conditions  with  the  efficiency  predicted  by 
the  model.  Unfortunately,  such  a  comparison  Is  beyond  the 
scope  of  this  study  (see  Appendix  B) • 

The  Problem 

It  Is  desired  to  design  a  Laser  oscillator  consisting 
of  a  rod  of  glass  doped  with  trlvalent  neodymium  and  pumped 
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by  an  annular  xenon  flash  lamp.  The  inner  surface  of  the 
lamp  Is  to  be  In  contact  with  the  rod,  and  under  these 
conditions  the  maxlmxim  energy  per  pulse  allowed  will  be 
given  by 

Eimax  =  Joules  (6-1) 

The  lamp  efficiency  in  the  range  .3  to  .9  microns  is  25^* 
and  the  geometrical  efficiency  of  the  pumping  scheme  is  12%, 
It  is  desired  to  maximize  the  energy  output  of  the  Laser, 
subject  to  the  following  restrictions: 

1 .  The  length  of  the  rod  will  be  25  cm. 

2.  The  minimum  radius  of  the  rod  is  0.3  cm. 

3.  The  maximxua  radius  of  the  rod  is  1.2  cm. 

4.  The  doping  density  is  7x10^^  lons/cm^. 

5.  The  reflectivity  of  the  output  end,  to  achieve  the 
desired  mode  selection,  must  be  at  least  50;^. 

Assume  that  the  ends  of  the  rod  are  to  be  coated  with 
multiple  layers  of  dielectric,  and  that  the  characteristics 
of  these  surfaces  are  given  by 

An  =  0.002Rn»  n  =  1,  2  (6-2) 

and  that  the  lamp  decay  time  constemt  'ty  is  given  by 

X]  =  AE^xIO”^  seconds  (6-3) 

Further  assxMe  that  the  lamp  spectrum  can  be  fairly  well 
represented  by  black  body  radiation  at  7000°K. 

A  sample  of  the  Laser  material  is  produced,  and 
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measurements  made  which  Indicate  the  fcllcwing  values  fcr 
the  physical  prcperties  cf  the  material: 

=  1.062  micrcns 
e  =  0.34 
idE  =  1930  cm“^ 

X  -  0.48  ms 
££.  =  1.19x10^^  cps 
0(0  =  0.015  nepers/cm 
cr(>0  as  given  in  lahle  I 


Table  I 

Abscrpticn  Crcss-secticns  cf  Katerial  in  10 


20 


cm 


X 

<r 

X 

C 

X 

(T 

X 

.30 

0. 

.46 

0.027 

.61 

0.300 

.76 

1 .460 

.31 

0. 

.47 

0.030 

.62 

0. 

.77 

0.352 

.32 

0. 

.48 

0.033 

.63 

0. 

.78 

0.278 

.33 

0. 

.49 

0. 

.64 

0. 

.79 

1.050 

.34 

0. 

.50 

0.200 

.65 

0. 

.80 

3.380 

.35 

0. 

.51 

0.740 

.66 

0. 

.81 

5.750 

.36 

0. 

.52 

0.860 

.67 

0.053 

.82 

1.300 

.37 

0. 

.53 

1 .810 

.68 

0.214 

.83 

0.374 

.38 

0. 

.54 

0.270 

.69 

0.181 

.84 

0.032 

.39 

0. 

.55 

0. 

.70 

0. 

.85 

0.032 

.40 

0. 

.56 

0.064 

.71 

0. 

.86 

0.342 

.41 

0. 

.57 

1 .300 

.72 

0.032 

.87 

0.779 

.42 

0. 

.58 

5.820 

.73 

0.370 

.88 

1.410 

.43 

0.117 

.59 

8.000 

.74 

3.230 

.89 

0.565 

.44 

0. 

.60 

1.880 

.75 

3.760 

.90 

0.235 

.45 

0.075 

The  values  cf  abscrptlcn  cross-section  In  Table  I  are 
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those  actually  measured  by  the  author  for  a  neodymlxun  doped 
lead  glass  sample  furnished  by  the  Isomet  Corporation, 
using  a  Cary  Model  14  spectrophotometer.  The  values  of  't 
and  /df  were  obtained  from  a  paper  describing  experiments 
with  neodymium  doped  barium  crown  glass  at  the  American 
Optical  Company  (Ref  2:364,  366),  The  value  of 
corresponds  to  a  fifty  percent  loss  In  about  46  cm.  The 
quantum  efficiency  Is  admittedly  a  guess  on  the  part  of  the 
author. 

If  It  Is  now  assumed  that  the  Laser  rod  will  be  operated 
at  room  temperature  (300°£),  enough  Information  Is  available 
to  compute  efficiency  as  a  function  of  radius  and  output  end 
reflectivity.  Multiplication  of  the  resulting  efficiencies 
by  the  appropriate  lamp  energies  will  give  the  output 
energies. 

Solution 

The  parameters  of  the  system  described  above  were 
entered  Into  the  computer  model  and  the  efficiency  and 
output  energy  calculated  for  ten  values  of  R  between  .3  cm 
and  1.2  cm,  and  for  five  values  of  R2  between  0.5  and  0.9* 

The  resulting  efficiencies  and  output  energies  are  listed  in 
Tables  II  and  III  on  the  following  page.  The  total  time 
required  to  compute  the  fifty  efficiencies,  including  the 
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time  required  to  warm  up  the  computer  and  load  the  programs, 
was  slightly  less  than  two  and  one-half  hours. 


Table  II 

Percent  Efficiency  vs.  Radius  and  Reflectivity. 


R(cm) 

R2=«5 

.6 

.7 

.8 

.9 

.5 

.286 

.241 

.193 

.137 

.073 

.4 

.335 

.282 

.224 

.159 

.084 

.5 

.373 

.315 

.251 

.178 

.095 

.6 

.405 

.341 

.273 

.193 

.103 

.7 

.428 

.363 

.290 

.206 

.110 

.8 

.447 

.382 

.306 

.218 

.116 

.9 

.462 

.396 

.317 

.226 

.121 

1.0 

“  “.765 

.328 

.234 

.125 

1.1 

1.2 

.446 

.708  \ 
.400 

-•236  _ 
.339  \ 

.242 

.246 

.129 

.133 

^Threshold  not  reached. 


Table  III 

Output  Energy  in  Joules  vs.  Radius  and  Reflectivity. 


R(cm) 

II 

Cn 

.6 

.7 

.8 

.9 

.3 

2.14 

1.81 

1.45 

1.03 

0.55 

.4 

3.35 

2.82 

2.24 

1.59 

0.84 

.5 

4.66 

3.94 

3.14 

2.23 

1.19 

.6 

6.08 

5.12 

4.09 

2.90 

.1.54 

.7 

7.49 

6.35 

5.08 

3.60 

1.93 

•  8 

9.94 

7.64 

6.12 

4.36 

2.32 

.9 

10.86 

9.31 

7.45 

5.31 

2.84 

1.0 

"  Tl76? 

\  10.18 

3.20 

5.85 

3.13 

1.1 

12.26 

“  ri72? 

6.65 

3.55 

1.2 

0.00 

12.00 

10.17  \ 

7.38 

3.99 
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An  examination  of  the  tables  shows  that  the  maximum  output 
energy  Is  obtained  when  the  radius  is  1.1  cm  and  the  output 
end  reflectivity  Is  0.5.  A  Laser  designed  using  these 
values  of  H  and  R2  should  produce  12.26  Joules  of  output 
with  an  overall  system  efficiency  of  0.446  percent.  Ihe 
maximum  efficiency,  on  the  other  hand,  Is  given  by  making 
the  radius  1.0  cm  and  the  output  end  reflectivity  0.5.  A 
Laser  built  using  these  values  should  produce  11.63  Joules 
of  output  with  an  efficiency  of  0.465  percent. 

One  feature  of  the  computer  model  which  Is  not  evident 
from  Tables  II  and  III  Is  that  the  computer  prints  the  value 
of  ,  the  time  to  threshold,  for  each  set  of  variables 
entered.  This  gives  the  operator  a  qualitative  Idea  of  how 
close  to  threshold  the  Laser  system  would  operate  under 
given  conditions.  In  the  sample  problem  above,  the  designer 
might  wish  to  maintain  a  safety  factor  to  insure  that  the 
Laser  constructed  would  actually  reach  threshold  even  though 
some  errors  might  be  present  In  the  measurement  of  the 
parameters  or  In  the  assumptions  made. 

One  way  to  do  this  Is  to  restrict  the  acceptable 
combinations  of  H  and  R2  to  those  which  reach  threshold  In 
less  than  a  certain  specified  time  period.  If,  as  an 
example,  no  solutions  are  considered  for  >riilch  T^  Is  greater 
than  200  microseconds,  the  combinations  of  R  and  R2  located 
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below  and  to  the  left  of  the  dashed  lines  In  Tables  II  and 
III  must  be  eliminated  from  consideration.  If  this 
restriction  Is  adopted,  the  maxlmxim  efficiency  (0.462^)  and 
maximum  output  energy  (10,86  J)  both  occur  when  H=0,9  and 
R2=0*5. 

In  the  sample  problem  above,  the  efficiency  and  output 
were  calculated  for  all  possible  combinations  of  the  two 
variables  being  optimized.  This  was  done  merely  to 
Illustrate  the  variation  In  efficiency  and  output  over  the 
entire  range  of  possible  combinations  of  R  and  R2*  In  sm 
actual  design  application.  It  would  probably  be  qulcher  and 
more  satisfactory  to  use  a  ''hlll-cllmblng"  method  of 
optimization,  particularly  if  several  parameters  were  to  be 
allowed  to  vary.  In  the  problem  above,  the  optimum  values 
of  R  and  R2  could  have  been  obtained  In  fifteen  tries  rather 
than  fifty,  even  If  the  worst  possible  starting  point  had 
been  selected  and  all  possible  mistakes  In  judgement  made. 
Care  must  be  exercised  In  the  use  of  this  method,  however, 
to  avoid  ending  on  top  of  the  wrong  "hill".  If  the  200 
microsecond  maximum  Is  Imposed  on  T^ ,  It  Is  possible  to 
obtain,  by  the  hlll-cllmblng  method,  an  optimum  solution  of 
R=1.2  and  R2=.8,  for  an  output  of  only  7.38  joules  at  an 
efficiency  of  0.246  percent.  Even  without  the  limit  on  Tp 
an  optlmxim  solution  of  R=1.2,  R2=.6  might  be  reached,  giving 
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an  output  of  12  Joules  at  0.4  percent  efficiency.  The  more 
closely  spaced  are  the  test  values  of  the  variables,  the 
less  is  the  chance  of  arriving  at  a  false  maximum  by  the 
hill-climbing  method. 

Summa-ry 

In  this  chapter,  the  use  of  the  computer  model  was 
Illustrated  by  the  solution  of  a  hypothetical  design  problem 
in  which  some  of  the  independent  variables  were  assigned 
fixed  values,  others  were  constrained  to  lie  within  certain 
limits,  and  still  others  were  made  functions  of  other 
variables.  The  flexibility  of  the  model  was  demonstrated 
by  optimizing  the  variables  with  respect  to  output  energy 
rather  than  efficiency,  and  by  introducing  the  possibility 
of  maintaining  a  safety  factor  to  insure  the  achievement  of 
threshold.  Fifty  trial  solutions  were  obtained,  and  the 
total  computer  operating  time  for  part  2  and  part  3  of  the 
program  was  less  than  two  and  one->half  hours.  The  use  of 
"hill-climbing"  techniques  in  the  optimization  process  was 
suggested  as  a  means  for  making  the  solution  even  faster, 
although  it  was  pointed  out  that  this  method  may  lead  to 
false  solutions  if  care  is  not  exercised. 

The  computer  printout  obtained  in  the  solution  of  the 
sample  problem  is  Included  in  this  paper  as  Appendix  C. 
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VII.  ReflultB  and  OonclueionB 

The  main  result  of  the  study  described  In  this  paper 
Is,  of  course,  the  mathematical  model  and  associated 
computer  program  for  the  determination  of  four-level  Laser 
efflolency.  Unfortunately,  time  limitations  have  prevented 
the  applloatlon  of  the  model  to  any  actual  problems  existing 
today  In  the  Laser  design  field.  Nonetheless,  some 
conclusions  may  be  dravn  from  the  model  and  from  the 
experience  gained  In  Its  develcpment  and  testing  regarding 
the  probable  effect  of  changing  each  of  the  independent 
variables  In  a  Laser  system  on  the  system  efflolenoy.  Some 
of  these  effects  have  already  been  physically  observed,  and 
this  fact  lends  support  to  the  validity  of  the  model.  In 
general,  it  oan  be  said  that,  within  the  limitations  imposed 
by  the  assiuaptions  and  approximations  used  in  the  develop¬ 
ment  of  the  model,  the  system  efficiency  will: 

1.  Increase  as  the  Laser  rod  length  is  increased. 

2.  Increase  as  the  Laser  rod  radius  is  increased. 

3*  Increase  as  the  doping  dexisity  is  increased. 

4.  Decrease  as  the  end  reflectivities  are  increased. 

5.  Increase  as  the  piusping  energy  is  increased. 

6.  Decrease  as  the  pumping  pulse  decay  time  constant 
is  increased. 
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Ho  attempt  has  been  made  to  make  the  model  applicable 
to  continuous  operation  of  the  four-level  Laser,  or  to  the 
so-called  pulsed  reflector  and  hair-trigger  modes  of 
operation.  It  Is  believed,  however,  that  the  basic 
mathematical  development  Is  applicable  to  all  of  these 
conditions,  and  that  only  minor  changes  In  the  computer 
program  would  be  necessary  In  order  to  be  able  to  represent 
them  by  the  model. 

Some  thought  was  given  to  the  possibility  of  making  the 
program  self-optimizing;  that  Is,  writing  the  program  In 
such  a  way  that  It  would  perform  automatically  the  hlll- 
cllmblng  optimization  referred  to  earlier  with  respect  to 
any  or  all  of  the  Independent  variables,  within  limits  set 
for  each  variable,  and  print  as  output  not  only  the  maximum 
obtainable  efficiency,  but  also  the  optimum  values  for  the 
variables.  It  was  felt,  however,  that  this  approach  could 
not  be  justified,  for  the  following  reasons: 

1 .  The  time  required  for  the  1 620  computer  to 
accomplish  an  adequate  optimization  by  the  hlll-cllmblng 
method,  particularly  If  many  variables  were  allowed  to  change 
and  small  Increments  of  each  variable  were  used,  would 
probably  be  quite  long. 

2.  It  would  be  very  difficult  to  make  the  program 
general  enough  to  allow  for  Interrelations  among  the 
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variables,  such  as  those  postulated  between  rod  dimensions 
and  lamp  parameters  In  the  sample  problem  of  Chapter  VI. 

3.  The  previously  mentioned  possibility  of  arriving 
at  a  false  optimum  point  through  the  use  of  hlll-ollmblng 
techniques  would  require  the  results  of  a  self-optimizing 
program  to  be  Interpreted  with  extreme  skepticism. 

4.  The  present  program  (part  3)  uses  most  of  the 
memory  of  the  1620,  and  It  Is  doubtful  that  self-optimization 
could  be  Included  In  a  program  lAilch  would  fit  In  a  20,000 
digit  memory. 

In  conclusion.  It  should  be  stressed  again  that  the 
computer  model  described  In  this  paper  has  not  been  oheoked 
against  any  experimental  results,  although  It  Is  known  that 
It  does  predict  the  correct  order  of  magnitude  of  efficiency 
and  threshold  for  neodymlvim  doped  glass  Lasers.  Before  the 
program  Is  used  In  actual  design  problems.  It  Is  strongly 
recommended  that  at  least  some  experimental  verification  of 
Its  adequacy  be  obtained.  This  project,  along  with  others 
which  the  author  feels  are  necessary  In  the  Laser  field.  Is 
Included  In  Appendix  B,  a  list  of  suggested  topics  for 
future  research. 
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Appaadix  A 


FORTRAN  Program 


PART  1 

PART  I,  4-LEVEL  LASER  EFFICIENCY  AND  GAIN  PR9GRAH  (IBH  1620). 
THIS  PROGRAM  CALCULATES  THE  CAPTIME  EFFICIENCY  VS.  RADIUS  TIMES 
ATTENUATION  CONSTANT. 

RUNNING  TIME  APPROXIMATELY  780  MINUTES 
DIMENSION  Fr20n,FRP(2n 
1  F0RMAT(E14.8) 

SR-SQRT(2.) 

BETA-O. 

00  101  1-1,201 

Ff I )-,015707963*EXPf-BETA) 

GAMMAa.031415927 
00  102  J«1,49 

F(l  )-Fr  I  )<t>.031 41  S927*EXPf -BETA/COS' GAMMA)) 

102  GAMMA-GAMMA4-.031415927 
101  BETAfBETA4.05 

ALPHR-.01 
00  110  J-1,21 
RP-C. 

00  109  K-1,21 

PHI-.01S707963 

FRP(K)-0. 

N-2 

A-.05 

00  106  L-1,100 

APR-SQRTf 1 .♦RP**2-2.*RP*C0S{ PHI ) ) 

BETA-ALPHR*APR 
00  124  l-N,201 
IF<BETA-A)12S,126,124 

1 25  FOFA-FI I )-( fi I )-F' 1-1 ) )* ( A-BETA)/.05 
GO  TO  127 

126  FOFAfFU) 

GO  TO  127 

124  A-A+.05 
GO  TO  104 

1 27  FRPf  X)-FRPf K )+.03l4l 5927*F0FA/APR 
1^1 

106  PHI-PHI<i>.03141S927 
104  FRP(K)-FRP(K)*ALPHR*RP 

109  RP^RP4,05 
MRELMCFRPr  1  )4>FRP  (21 )  )/2. 

DO  103  1^2,20 

103  NREL/^WRELAMTRPCN) 

MRELA-WRELA/31 .41 5927 
PUNCH  I.MRELA 

110  ALPHN-ALPHR*SR 
STOP 

END 
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PART  2 

PART  2,  4-LEVEL  laser  EFFICIENCY  AND  GAIN  PR9GRAM  (IBM  1620). 
THIS  PROGRAM  CALCULATES  THE  PRODUCT  OF  CAPTURE  EFFICIENCY  AND 
RELAXATION  EFFICIENCY  VS.  RADIUS  TIMES  DOPING  DENSITY, 

RUNNING  TIME  APPROXIMATELY  25  MINUTES. 

DIMENSION  Hi:i01),SIGMA(l01),F(101),Wf  10n,WRELA(2l) 

1  F0RMAT:e14.8) 

2  FORMAT (/41HSWITCH  I  ON  TO  REPRESENT  LAMP  SPECTRUM  BY) 

3  FORMAT' 3411  BLACK  BODY  RADIATION,  PUSH  START./) 

4  FORMAT :/42HENTER  BLACK  BODY  TEMPERATURE  IN  DEGREES  K.) 

5  FORMAT ;/47HENTER  LIMITS  OF  ABSORPTION  SPECTRUM  IN  MICRONS.) 

6  format ;/4IHENTER  laser  OUTPUT  WAVELENGTH  IN  MICRONS.) 

7  FORMAT ;/25HENTER  QUANTUM  EFFICIENCY.) 

PRINT  2 

PRINT  3 
PAUSE 
PRINT  7 
ACCEPT-E 
PRINT  5 

ACCEPT, WV1,WV2 
IF,WV1-WV2)101.101,102 

102  X-WV1 
WVt-WV2 
WV2-X 

101  PRINT  6 
ACCEPT,  W*/ 

X-WV1 

IFiSENSE  SWITCH  1)103,104 

103  PRINT  4 
ACCEPT, T 

DO  105  1-1,101 
READ, sigma' I ) 

Hf|)-(X**(-5))/(EXPM4496.35/a*X))-1,) 

X-X+.OI  . 

IF(X-WV2)10S,105,106 

105  CONTINUE 

104  DO  107  1-1,101 
READ,SIGMA(I),H(I) 

X-X-I-.01 

IF(X-WV2)107,107,106 

107  CONTINUE 

106  SUM-^H( l)+H(1))/2. 

J-l-1 

DO  108  K-2,J 

108  SUFWiSUM4-H'K) 

00  109  K-1,1 

109  H(K)-H(K)*'WV2-WV1)/^01*SUH) 

00  110  K-1,21 

110  READ  I.WRELAVK) 

SR-SQRT(  2. ) 

RH0R-.OO1 

DO  111  K-1,21 
00  112  L-1,1 
112  Fa)-RH0R*SIGMA(.L) 

X-WV1 

00  120  L-1,1 
ALPHR-.01 
00  113  H-t,2l 
IFIALPHR-FrL))113, 114.115 
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113  ALPHR-ALPHR*SR 

114  W(L)-WRELA<M)*Ha)*X 
Gft  T«  120 

115  IF(M-nil7,117,1l8 

117  WfL)-0, 

Gt  T0  120 

118  WfL)-WRELAfM)-a0G(ALPHR/F(L))*(WRELA(M)-WRELAfM-n)/LM(SR)) 

wa)-wa)*HfL)*x 


120 


119 


111 


X-X4-,01 

SUlVi'wrn+WM))/2. 

Of  119  L-2,J 
SlAViSUKfWa) 

ETA-.Ol^SUM^E/ 'WV*(WV2-WV1 
PUNCH  l.ETA 
RH0R-RH0R*SR 


)) 


STOP 

END 


PART  3 

PART  3.  4>LEVEL  LASER  EFFICIENCY  AND  GAIN  PROGRAM  (IBM  1620). 

THIS  ptoAM  CALCULATES  THE  TIME  Ti  THRESHBLD  AND  EFFICIENCY  BF  A 
LASER  OSCILLATOR  ANO  THE  TIME  TO  THRESHOLD, MAXIMUM  EFFICIENCY  AND 
MAXIMUM  GAIN  OF  A  LASER  AMPLIFIER. 

RUNNING  TIMES  APPROXIMATE)  >  PRECOMPUTATION.  2  MINUTES. 
OSCILLATOR  PROBLEM,  3  MINUTES.  AMPLIFIER  PROBLEM.  2  MINUTES. 
DIMENSION  ETA(21).GAMMA(32) 

1  FORMAT :/34HENTER  ROD  RADIUS  AND  LENGTH  IN  CM.) 

F0RMAT(/33HEMTER  REFLECTIVITY  OF  OUTPUT  END.) 

F0RMAT(/51HENTER  lamp  energy  in  joules  ANO  DECAY  TIME  CONSTANT) 
F0RMAT'17H  in  MICROSECONDS.) 

5  FORMAT (/36HENTER  DOPING  DENSITY  IN  IONS  PER  CC.) 

FORMAT ;/41HEMTER  OPERATING  TEMPERATURE  IN  DEGREES  K.) 

FORMAT ;/45HSWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR) 
F0RMATn3H,  PUSH  START./) 

FORMAT (/45HEMTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END.) 

10  FORMAT A/ 19HRH0  R  OUT  OF  RANGE./) 

11  FORMAT(//44HTHRESHOLO  NOT  REACHED,  INCREASE  LAMP  ENERGY,/) 

12  FORMAT <;/44HEHTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL.) 

|j  F0RMAT7/19HTIME  TO  THRESHOLD  -,F5.0,14H  MICROSECONDS,) 

II  max.  gain -.F12.7) 

7  F0^r/46lCNTER  LAMP  EFFICiInCY.  iEOMETRICAL  bFICIENCY.) 

18  F0^T(/5OHENTER  FLUORESCENT  LInI  WIDTH  IN  CVaES  PER  SECOND.) 

9  FORMATf/llHEMTER  LIFETIME  OF  METASTABLE  STATE  IN  MICROSECONDS.) 

20  FOWATl/MHENTER  ENERGY  DIFFERENCE  OETVCEN  GROUND  STATE  AND) 

21  FORMAT ;33H  TERMINAL  STATE  IN  RECIPROCAL  CM.) 

22  FORMAT (/41HENTER  LASER  OUTPUT  WAVELENGTH  IN  MICRONS.) 

23  format;///) 

DO  101  1^1,21 

101  READ  13.ETa;n) 

print  20 

PRINT  21 
ACCEPT.DEL 
PRINT  22 
ACCIPT,WV 
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PRINT  19 
ACCEPT, TAU 
PRINT  17 

ACCEPT, ETAL.ETAG 

A— ETAL*ETAG*WV*TAU/19.86 

T-20. 

D«  100  1^1,32 
GAMMAOO-EXPUT/TAU) 

100  T-T+20. 

PRINT  10 
accept, OELF 
PRINT  6 
ACCEPT, TK 
PRINT  23 
123  PRINT  12 

ACCEPTjALPHA 
10S  PRINT  5 
ArcFPT  R>4a 

P0P-RHi*EXP'-1  .i»496*DEL/TK) 
PRINT  1 
accept, R,RL 
RH«R-1.E-2O*RH0*R 
IFfRHaR-.001)108.109,109 
109  IF(1.024-RHOR)108,104,104 
108  PRINT  10 
G«  T«  105 
104  PRINT  7 
PRINT  8 


PAUSE 

S I GMA- 3 . 74E-04* ( WV**2 ) / < OELF*TAU ) 

IF^SENSE  SWITCH  1)114,115 

114  PflP«P0P-»-ALPHA/SIGHA 
G0  T0  116 

115  PRINT  2 
ACCEPT, R2 
PRINT  9 
ACCEPT, A1,A2 

R1-1.-A1  .  , 

P0P-P0P-» ;  2.*ALPHA*RL-L0G  (R1^2))/(2.*SI  GMA*ftL ) 

116  PRINT  3 
PRINT  4 
ACCEPT, El, TAU1 
Z-.001*SQRTr2.) 

00  117  1^2,21 
IF^RH0R«Z)1I8,119,117 

117  Z-Z*SQRTf2,) 

119  E-ETA'N) 

G0  T0  120 

1 18  E-ETAf  lO-LOCr  Z/RH0R)*' ETA^  N)-ETA' N-1 ) )/L0Gr  SQRT^  2. ) ) 

120  TAU>TAU1/10.35 
0EL1-TAU1-TAU2 
0EL2-TAU1/'TAU-TAU1) 

0EL3-TAU2/'TAU-TAU2) 

DELV0EL24>DEL3 

T-20. 

00  121  1^1,32 

P3-A*rOEL2*EXPUT/TAU1)-i>DEL3*EXP(-T/TAU2)-DEL4*6AMHA(N))/OEL1 
P3-1  .E+20*P3*E1*E/  •  3. 1415927*'A**2)*IU.) 

IFfP3-P8P)t21.122,l22 
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121 

122 

103 

102 

129 

130 

106 


T-T+20, 

PRINT  11 
G«  T0  116 
PRINT  14, T 

EFF-TAU1*0EL2*EXP (-T/TAU1 )+TAU2*0EL3*EXP(-T/TAU2) 

EFF— 1.E-»02*(EFF-TAU*OEL4*GAHHA'N))/DEL1*ETAL*ETAG*E 
IFISENSE  SWITCH  1)102,103 
GNET— L06fR1*R2)/f2.*RL) 

0EN-GHET*M,-R2) 

0EI^DEN*'ALPHA*R2*<:  1 .4R1*EXP' 6NET*RL)  )*' EXPr  6NET*RL)-1  •  ) 
DEN-0EIHA1*R2*GNET*EXP^GNET*RL) 

EFF-EFFfGNET*M.-R2-A2)/0EN 

PRINT  fS.EFF 

G0  T0  106 

0«  129  K-N,32 

F-P3 

P>A*r  0EL2*EXPr-T/TAU1  )«0EL3*|XP'-T/TAU2)-0EL4*GAMHA^IC)  )/0EL1 
P> 1 . E+20*P3*E 1*E/ ^3.141 5927*(R**2)*RL ) 

IFrP3-Fl130,129.129 

T-T+20. 

F-P3 

G-EXP<RL*SIGMA*fF>P0P)) 

EFF-EFF/< 1 .+ALPHA*RL) 

PRINT  16,EFF,G 
PRINT  23 
PAUSE 
G0  T0  123 
tNO 
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Appendix  B 


Suggested  Topics  for  Future  Research 


1 .  The  development  of  accurate  experimental  methods  for  the 
measurement  of  the  physical  properties  of  Laser 
materials . 

2.  The  measurement  and  tabulation  of  the  physical 
properties  of  Laser  materials  under  various  conditions. 

3.  An  investigation  of  the  efficiency  of  present  pumping 
sources  and  Laser  pumping  geometries. 

4.  A  detailed  study  of  the  spectral  output  of  Laser  piunplng 
sources. 

5.  An  experimental  check  of  the  results  of  the  computer 
model  developed  In  this  paper. 

6.  An  extension  of  the  applicability  of  the  computer  model 
to  elliptical  pximplng  geometries,  and  to  the  continuous, 
pulsed  reflector,  and  hair-trigger  modes  of  operation. 


71 


GE/EE/62-18 


Appendix  2, 


Computer  Printout  for  Sample  Problem 


ENTER  ENERGY  DIFFERENCE  BETWEEN  GROUND  STATE  AND  TERMINAL  STATE  IN  RECIPROCAL  CM. 
1930 

ENTER  LASER  OUTPUT  WAVELENGTH  IN  MICRONS. 

1.062 

ENTER  LIFETIME  OF  METASTABLE  STATE  IN  MICROSECONDS. 

480 

ENTER  LAMP  EFFICIENCY.  GEOMETRICAL  EFFICIENCY. 

.25  .72 

ENTER  FLUORESCENT  LINE  WIDTH  IN  CYCLES  PER  SECOND. 

1.19E13 

ENTER  OPERATING  TEMPERATURE  IN  DEGREES  K. 

300 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.3  25 

SWITCH  1  ON  FOR  AMPLIFIER.  OFF  FOR  OSCILLATOR.  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END.  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS^ 
750  300 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS. 

EFFICIENCY  -  .28623733  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.3  25 

SWITCH  1  ON  FOR  AMPLIFIER.  OFF  FOR  OSCILLATOR.  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  absorptivity  OF  CLOSED  END.  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
750  300 

TIME  TO  THRESHOLD  -  UO.  MICROSECONDS. 

EFFICIENCY  -  .24115724  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.3  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 U 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
750  300 

TIME  TO  THRESHOLD  -  20.  MICROSECONDS, 

EFFICIENCY  -  .19256235  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.3  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  EMU,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
750  300 

TIME  TO  THRESHOLD  -  20.  MICROSECONDS, 

EFFICIENCY  -  .13652181  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.3  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
750  300 

TIME  TO  THRESHOLD  -  20.  MICROSECONDS, 

EFFICIENCY  -  .0725765'*  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.'t  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1000  400 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS. 

EFFICIENCY  -  .33482683  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PERCC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.4  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1000  400 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS. 

EFFICIENCY  -  .28209428  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.4  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.7 

ENTER  absorptivity  OF  CLOSED  END,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1000  400 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS, 

EFFICIENCY  -  .22397586  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.4  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.8 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1000  400 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS, 

EFFICIENCY  -  .15877670  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.4  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1000  400 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS. 

EFFICIENCY  -  .08440750  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC, 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

,5  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START, 

ENTER  reflectivity  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END, 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1250  500 

TIME  TO  THRESHOLD  -  60.  MICROSECONDS, 

EFFICIENCY  -  .37345961  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.5  25 

S\nTCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1250  500 

TIME  TO  THRESHOLD  -  60.  MICROSECONDS, 

EFFICIENCY  -  .31464271  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROO  RADIUS  AND  LENGTH  IN  CM. 

.5  25 

SWITCH  1  OH  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1250  500 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS, 

EFFICIENCY  -  .25124lt05  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

n\'d 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.5  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.8 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1250  500 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS, 

EFFICIENCY  -  .17810714  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E1‘J 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.5  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1250  500 

TIME  TO  THRESHOLD  -  AO.  MICROSECONDS, 

EFFICIENCY  -  .0‘JA68T77  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AMD  LENGTH  IN  CM. 

.6  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1500  600 

TIME  TO  THRESHOLD  -  80.  MICROSECONDS, 

EFFICIENCY  -  .40501937  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.6  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END, 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1500  600 

TIME  TO  THRESHOLO  -  00.  MICROSECONDS, 

EFFICIENCY  -  .34123207  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.6  25 

r,\<ITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1500  600 

TIME  TO  THRESHOLD  -  60.  MICROSECONDS, 

EFFICIENCY  -  .27264604  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC; 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.6  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.fl 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1500  600 

TIME  TO  THRESHOLD  -  60.  MICROSECONDS, 

EFFICIENCY  -  .19T27961  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.6  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START, 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1500  600 

TIME  TO  THRESHOLD  -  40.  MICROSECONDS, 

EFFICIENCY  -  .10325479  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.7  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.OD2  .001. 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS 
1750  700 

TIME  TO  THRESHOLD  -  120.  MICROSECONDS, 

EFFICIENCY  »  .4271)6621  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.7  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  ANO  DECAY  TIME  CONSTANT  IN  MICROSECONDS 
1750  700 

TIME  TO  THRESHOLD  -  100.  MICROSECONDS, 

EFFICIENCY  -  .363^^719  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.7  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  ANO  DECAY  TIME  CONSTANT  IN  MICROSECONDS 
1750  700 

TIME  TO  THRESHOLD  -  80.  MICROSECONDS, 

EFFICIENCY  -  .29041761  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  miST  MATERIAL. 
.015 

ENTER  doping  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.7  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START, 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.0 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1750  700 

TIME  TO  THRESHOLD  -  80.  MICROSECONDS, 

EFFICIENCY  -  .205C7731  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E13 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.7  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  ENO,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1750  700 

TIME  TO  THRESHOLD  -  60.  MICROSECONDS, 

EFFICIENCY  -  .11006010  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.0  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  reflectivity  OF  OUTPUT  ENO. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2000  800 

TIME  TO  THRESHOLD  -  160.  MICROSECONDS. 

EFFICIENCY  -  .44654105  PERCENT. 


ENTER  ATTENUATIO>l  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7EI9 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 


80 


Gii/3ii/62-1  3 


8  23 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 

2000  Boo 

TIME  TO  THRESHOLD  -  120.  MICROSECONDS, 

EFFICIENCY  -  .3B233702  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.8  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 

2000  800 

TIME  TO  THRESHOLD  -  100.  MICROSECONDS, 

EFFICIENCY  -  .30567409  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

8  23 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START, 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.0 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2000  ROO 

TIME  TO  THRESHOLD  ■  OO.  MICROSECONDS, 

EFFICIENCY  .  .21800274  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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.n  25 

S'/ITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START, 

r:.:tk;i  reflectivity  of  output  end. 

EilTER  ARC  RRTIVITY  OF  CLOSEO  END,  OUTPUT  END. 

.002  .001 r 

EIITER  LAMP  EllERCY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS, 
2000  f.OO 

TIME  T!'  THRESH'OLO  -  CO.  MICROSECONDS, 

EFFICIEIICY  -  .115r.y273  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 

.9  25 

SWITCH  1  ON  FOR  AMPLIFIER.  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END, 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  ANO  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2250  900 

TIME  TO  THRESHOLD  -  200.  MICROSECONDS, 

EFFICIENCY  -  ,461753*»7  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 

.9  25 

S'./ ITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START, 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2250  900 

TIME  TO  THRESHOLD  -  160.  MICROSECONDS, 

EFFICIENCY  -  .39573184  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7F19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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/;  25 

r.'JITW  1  UN  Ft;R  AMPLIFIER,  (IFF  F('.R  LSCI LLAT(:R,  PUSH  START, 

EilTF.R  REFLECTIVITY  (.F  {'UTPIJT  Elll). 

.7 

EllTER  ASS'/RPTIVITY  SF  CL('SEl)  EHi),  OUTPUT  EMU. 

.002  .ooin 

EUTER  LAMP  EMERSY  111  .KHILES  AMD  DECAY  TIME  CCMSTAHT  IN  Ml  CRdSECf'HDS 

2250  5:00 

TIME  Tf‘.  THRESHCLl)  »  lUO.  Ml CROSECIitlOS, 

EFFICIENCY  -  .3l66771't  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  OCPIMG  DENSITY  IN  IONS  PER  CC. 

7C1‘J 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

.‘J  25 

SWITCH  1  (!N  FOR  AMPLIFIER,  (‘FF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  t'F  CllTPlIT  Elll). 

.R 

ENTER  ADSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERCY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS 
2250  DOO 

TIME  TO  THRESHOLD  -  120.  MICROSECONDS, 

EFFICIENCY  -  .226o62M(  PERCENT. 


ENTER  ATTENUATION  C!  NSTANT  ('F  HOST  MATERIAL. 

.015 

ENTER  OOP  I  NO  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 

.9  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F(’R  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ADSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  ANO  DECAY  TIME  CONSTANT  IN  MICROSECONDS 
2250  900 

TIME  TO  THRESHOLD  -  100.  MICROSECONDS, 

EFFICIENCY  -  .1209'i952  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 
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1  25 

SWITCH  I  0N  FOR  AMPLIFIER.  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.5 

ENTER  ASSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2500  1000 

TIME  TO  THRESHOLD  -  2fl0.  MICROSECONDS, 

EFFICIENCY  -  .46515602  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1  25 

SWITCH  I  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2500  1000 

TIME  TO  THRESHOLD  «  200.  MICROSECONDS. 

EFFICIENCY  -  .40678021  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  ENO,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2500  1000 

TIME  TO  THRESHOLD  -  160.  MICROSECONDS, 

EFFICIENCY  -  .H2815479  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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1  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START. 

F.NTER  REFLECTIVITY  OF  OUTPUT  END. 

.8 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2500  1000 

TIME  TO  THRESHOLD  -  140.  MICROSECONDS, 

EFFICIENCY  -  .23429432  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 

1  25 

SWITCH  1  ON  FOR  AMPLIFIER.  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  ANO  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2500  1000 

TIME  TO  THRESHOLD  -  120.  MICROSECONDS, 

EFFICIENCY  -  .12536659  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.1  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START, 

ENTER  REFLECTIVITY  OF  OUTPUT  END, 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2750  1100 

TIME  TO  THRESHOLD  -  440.  MICROSECONDS, 

EFFICIENCY  -  .44557654  PERCENT, 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 
.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 
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1.1  25 

SWITCH  1  'IN  FOR  AMPLIFIER,  0FF  FAR  0SCILLAT0R,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  ENO. 

.6 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  ENO. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2750  1100 

TIME  TO  THRESHOLD  -  280.  MICROSECONDS, 

EFFICIENCY  -  .4081(7202  PERCENT. 


ENTER  ATTENUATION  CONSTANT  0F  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.1  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0014 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2750  1100 

TIME  TO  THRESHOLD  -  200.  MICROSECONDS, 

EFFICIENCY  «  .33578956  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.1  25 

SWITCH  1  ON  FuR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  ENO. 

.8 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  ENO. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2750  1100 

TIME  TO  THRESHOLD  -  160.  MICROSECONDS, 

EFFICIENCY  -  .24159533  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  ANO  LENGTH  IN  CM. 
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1.1  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  F0R  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
2750  1100 

TIME  TO  THRESHOLD  -  140.  MICROSECONDS, 

EFFICIENCY  -  .12928915  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.'2  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.5 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.<002  .001 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
3000  1200 

THRESHOLD  NOT  REACHED,  INCREASE  LAMP  ENERGY. 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
1E6  1000 

TIME  TO  THRESHOLD  -  20.  MICROSECONDS, 

EFFICIENCY  -  .54229826  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.2  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0012 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
3000  1200 

TIME  TO  THRESHOLD  -  400.  MICROSECONDS, 

EFFICIENCY  -  .39978959  PERCENT. 
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ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 


ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.?.  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 


ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.7 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.0D2  .001l» 

ENTER  LAI^  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
SOOO  1200 

TIME  TO  THRESHOLD  -  260.  MICROSECONDS, 

EFFICIENCY  -  .ST931019  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E19 

ENTER  Ri'J  RADIUS  AND  LENGTH  IN  CM. 

1.2  25 

SWITCH  1  1)11  F;jR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 


ENTER  REFLECTIVITY  OF  OUTPUT  END. 


ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0016 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
BOOO  1200 


TIME  TO  THRESHOLD  -  200.  MICROSECONDS, 
EFFICIENCY  -  .24636118  PERCENT. 


ENTER  ATTENUATION  CONSTANT  OF  HOST  MATERIAL. 

.015 

ENTER  DOPING  DENSITY  IN  IONS  PER  CC. 

7E1‘J 

ENTER  ROD  RADIUS  AND  LENGTH  IN  CM. 

1.2  25 

SWITCH  1  ON  FOR  AMPLIFIER,  OFF  FOR  OSCILLATOR,  PUSH  START. 

ENTER  REFLECTIVITY  OF  OUTPUT  END. 

.9 

ENTER  ABSORPTIVITY  OF  CLOSED  END,  OUTPUT  END. 

.002  .0018 

ENTER  LAMP  ENERGY  IN  JOULES  AND  DECAY  TIME  CONSTANT  IN  MICROSECONDS. 
3000  1200 

TIME  TO  THRESHOLD  -  160.  MICROSECONDS, 

EFFICIENCY  -  .13279791  PERCENT. 
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